94-29881 


^8 


February  1994 


Prepared  for 

Ballistic  Missile  Defense  Organization 


Approved  tor  public  nlutt;  dlttrltwtlon  onllmttod. 


//b-fS'M  /£ 

Copy  ]  2  of  45  copies 


AD-A284 

IVHH111I 


IDA  DOCUMENT  D-1483 


IDA  GAMMA-RAY  LASER 
1993  ANNUAL  REPORT 


Bohdan  Balko 
Irvin  W.  Kay 


DTIC 

ELECVitl  jp'ji 

St?  1  5  ,994  |S  !] 


DA 


iync  QUALITY  INSPECTED  3 

INSTITUTE  FOR  DEFENSE  ANALYSES 

1801  N.  Beauregard  Street,  Alexandria,  Virginia  22311-1772 


94  9  14  023 


IDA  Log  No.  HQ  94-45195 


Best 

Available 

Copy 


DEFINITIONS 

IDA  publishes  the  following  documents  to  report  too  recalls  of  Ns  work. 

Reports 

Reports  are  too  most  authoritative  anO  most  carefully  considered  products  IDA  publishes. 
They  normally  embody  results  ot  mo|or  prelects  which  (»)  hays  a  Street  bearing  on 
Seclslons  affecting  major  prof  rams,  (b)  aSSrass  lasaas  ol  significant  cancam  to  toe 
ExacoUva  Branch,  too  Congress  aaS/ar  toe  pabllc,  or  (c)  aSSrass  lasaas  that  have 
significant  economic  implications.  IDA  Reports  are  reviewed  by  ovtside  panels  of  experts 
to  enema  their  high  gaallty  anS  relevance  te  toe  problems  sbiSleS,  ana  they  are  releases 
bythePresISentetIDA. 

Group  Reports 

Group  Reports  recorf  the  flnSIngs  anS  resells  ol  IDA  establishes  working  groups  anS 
panels  composes  ot  senior  InShrtfoals  addressing  major  Issaea  which  otherwise  wools  he 
toe  cabled  el  an  IDA  Report.  IDA  Group  Reports  are  revtewef  by  toe  senior  inSMSeals 
responsible  lor  the  protect  anS  others  as  selectef  by  IDA  to  ensure  their  high  quality  anS 
relevance  to  the  problems  stuSleS,  and  are  releases  by  toe  President  ol  IDA. 


Papers 

Papers,  else  authoritative  anS  carefully  consif erof  profects  ol  IDA,  address  studies  that 
are  narrower  In  scope  than  those  covered  In  Reports.  IDA  Papeis  are  reviewed  to  ensure 
that  they  meet  toe  high  standards  expected  of  refereed  papers  in  professional  journals  or 
formal  Agency  reports. 


Documents 

IDA  Documents  are  used  for  the  convenience  ol  the  sponsors  or  toe  analysts  (a)  to  record 
substantive  work  done  In  quick  readlon  studies,  (b)  to  record  the  proceedings  of 
conferences  and  meetings,  (c)  to  make  available  preliminary  and  tentative  results  ol 
analyses,  (d)  la  record  data  developed  In  the  comae  ot  an  Investigation,  or  (a)  to  forward 
Information  that  Is  essentially  unanalyzed  and  unevaluated.  The  review  of  IDA  Documents 
is  suited  to  their  content  end  Intended  use. 


The  work  reported  In  tola  document  was  conducted  under  contract  MDA  903  89  C  0003  lor 
the  Department  of  Defence.  The  peblf cation  of  tola  IDA  document  does  not  indicate 
endorsement  by  the  Department  ot  Defense,  ner  should  the  contents  be  construed  as 
reflecting  toe  official  position  ot  that  Agency. 


IDA  DOCUMENT  D-1483 


IDA  GAMMA-RAY  LASER 
1993  ANNUAL  REPORT 


#  INSTITUTE  FOR  DEFENSE  ANALYSES 

Contract  MDA  903  89  C  0003 
Task  T-R2-597.03 


DTIC  QUALITY  1_L1S?ECTED  3 


0.0 


PREFACE 


t 


This  document  represents  work  done  by  the  Institute  for  Defense  Analyses  in 
FY  1993  under  Task  Order  T-R2- 597.3,  Technical  Support  for  Innovative  Concepts 
Program,  for  the  Ballistic  Missile  Defense  Organization. 


5 


ii 


ABSTRACT 


► 


I 


This  document  consists  of  a  collection  of  papers  published  in  the  open  literature 
1  in  FY  1993. 

The  first  paper  was  published  in  Phys.  Rev.  B  and  describes  the  limitations  of  the 
Bonifacio-Lugiato  model  of  superfluorescence. 

|  The  second  paper  deals  with  the  feasibility  of  nuclear  superfluorescence.  It  was 

presented  at  the  Lasers  *92  Conference  in  Houston  and  published  in  the  Proceedings. 

Paper  Number  3  describes  a  nuclear  superfiuorescence  model.  It  was  presented  at 
the  Physics  in  Ukraine  Conference  in  Kiev  in  June  1993  and  published  in  the  Proceedings. 

» 

Paper  Number  4  is  a  comment  on  a  paper  published  in  Phys.  Rev.  B,  which  we 
found  to  be  wrong  and  which  stated  a  misleading  conclusion  regarding  the  effect  of 
homogeneous  broadening  on  an  inhomogeneously  broadened  system.  Our  comment  is 
^  followed  by  the  author's  reply,  which  was  also  published  in  the  Comments  section  of 

Phys.  Rev.  B.  We  include  our  response  to  the  author's  reply. 
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Critique  of  the  Bonifacio-Lugiato  superfluorescence  model 

I.  W.  Kay  and  B.  Balko 

Institute  for  Defense  Analyses.  Alexandria.  Virginia  22311 
(Received  23  October  1992) 

This  paper  shows  bow  the  Bonifscio-Luguto  mean-field  model  of  atomic  superfluorescence,  which  is 
particularly  attractive  because  it  deals  with  a  crystalline  array  of  emitters,  can  be  extended  to  the  nuclear 
case.  However,  limits  exist  on  the  validity  of  certain  approximations  used  in  the  development  of  the 
model.  In  particular,  a  numerical  instability  that  increases  without  bound  is  shown  to  exist  in  the  non- 
Markovian  case,  placing  in  doubt  the  application  of  the  semiclaasical  theory  derived  for  the  non- 
Markovian  case. 


L  INTRODUCTION 

Recent  activity  in  the  field  of  y-ray  lasers  has  focused 
on  the  possibility  of  observing  nuclear  superfluorescence 
(SF).  The  most  general  treatment  of  SF  based  on  the 
Maxwell-Bloch  (MB)  equations'  permits  consideration  of 
several  important  phenomena2  which,  while  interesting 
only  as  second-order  effects  in  atomic  SF,  are  crucial  to 
the  understanding  of  nuclear  superfluorescence.  Con¬ 
sideration  of  emitters  in  a  crystal  lattice,  which  leads  to 
special  effects  such  as  Borrmann  anomalous  transmis¬ 
sion,2  are  not  treated  in  the  MB  formalism.  On  the  other 
hand  the  model  that  Bonifacio  and  Lugiato  presented  in 
1973  (Refs.  4  and  3)  deals  with  emitters  in  a  crystal  lattice 
and  may  be  adaptable  to  a  similar  structure  of  nuclear 
emitters  that  radiate  electromagnetic  energy  at  wave¬ 
lengths  smaller  than  the  interatomic  spacing.  This  possi¬ 
bility  and  the  quantum-mechanical  basis  of  the  BL  model 
motivated  us  to  revisit  the  early  work.  This  paper 
presents  some  of  our  findings  regarding  the  BL  model. 

The  BL  theoretical  model  for  the  behavior  of 
superfluorescence  by  a  system  of  identical  two-level 
atoms  remains  an  important  contribution  to  a  fundamen¬ 
tal  understanding  of  that  phenomenon,  despite  the  fact 
that  the  theory  does  not  properly  account  for  propaga¬ 
tion  or  the  initial  quantum  fluctuations  of  the  electric- 
field  vacuum  state.  For  example.  Polder,  Schurmans,  and 
v  rehen6  developed  a  model  of  a  similar  nature  without 
these  defects;  however,  its  derivation  requires  an  assump¬ 
tion  that  the  atoms  occupy  random  positions,  thereby 
ruling  out  its  application  to  a  crystal  lattice.  But  not  only 
does  the  BL  model  specifically  assume  that  the  radiating 
structure  is  a  crystal,  it  also  takes  into  account  inhomo- 
geneous  broadening  due  to  motion  of  the  atoms  about 
their  periodically  spaced  equilibrium  positions,  whereas 
Ref-  6  neglects  inhomogeneous  broadening. 

The  crystalline  structure  is  an  important  consideration 
when  nuclear  transitions  resulting  in  photon  wavelengths 
smaller  than  the  interatomic  spacing  are  responsible  for 
the  superfluorescent  emission.7  Although  in  deriving  the 
BL  model,  the  authors  assume  only  atomic  transitions  for 
which  the  photon  wavelengths  are  larger  than  the  intera¬ 
tomic  spacing,  the  associated  analysis  is  easily  adapted  to 
the  nuclear  case.  Section  II  will  show  that  when  the 
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model  is  applied  to  the  case  of  radiation  wavelengths 
smaller  than  the  interatomic  spacing  it  predicts  that 
superfluorescence  is  then  only  possible  at  wavelengths 
satisfying  certain  conditions.  It  turns  out  that  these  con¬ 
ditions  are  also  sufficient  for  coupling  to  a  Borrmann 
propagation  mode,  which  has  been  cited  as  necessary  for 
y-ray  lasing.1 

Although  the  BL  model  approaches  a  classical  limit  in 
the  initial  stage,  during  which  certain  characteristics  of 
the  emitted  pulse,  such  as  the  time  delay,  are  determined, 
the  treatment  is  quantum  mechanical.  Basing  calcula¬ 
tions  on  the  more  accurate  quantum-mechanical  equa¬ 
tions,  Sec.  Ill,  in  fact,  presents  numerical  results  demon¬ 
strating  that  the  classical  approximation  may  not  be  valid 
until  well  after  the  radioed  pulse  maximum  has  emerged. 
However,  Sec.  IV  shows  that  in  certain  circumstances  a 
numerical  instability  invalidating  the  quantum- 
mechanical  calculations  will  occur. 

II.  short-wavelength  superfluorescence 

It  wavelengths  associated  with  nuclear  transitions  are 
smaller  than  interatomic  spacings,  efficient  radiation  can 
only  occur  when  conditions  for  Borrmann  mode  propaga¬ 
tion  are  satisfied.  Since  superfluorescence  is  also  deemed 
necessary  for  lasing  in  the  y-ray  regime,  the  effect  of 
wavelengths  shorter  than  interatomic  spacings  on 
superfluorescence  is  of  interest.  In  particular,  it  is  of 
some  interest  to  determine  the  extent  to  which  conditions 
for  superfluorescence  at  short  wavelengths  may  be  in 
conflict  with  those  for  Borrmann  mode  propagation. 

In  nuclear  superfluorescence,  wavelengths  associated 
with  transitions  of  interest  will  undoubtedly  be  smaller 
than  the  spacing  of  the  contributing  nuclei  in  a  crystal 
lattice.  The  behavior  of  a  similarly  spaced  classical 
periodic  array  of  coherent  sources  suggests  that  for  a  po¬ 
tentially  superfluorescent  crystal,  regarded  as  an  array  of 
identical  nuclei,  a  first-order  effect  of  such  a  short  wave¬ 
length  would  be  to  limit  collective  radiation  to  certain 
directions.  In  fact,  the  relatively  simple  BL  model  of 
atomic  superfluorescence  has  sufficient  scope  to  confirm 
the  existence  of  such  an  effect. 

The  BL  superfluorescence  model  for  two-state  atoms  in 
a  rectangular  crystal  lattice  considers  only  the  case  in 
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which  the  emitted  radiation  wavelength  is  larger  than  the 
atomic  spacing.  However,  this  limitation  is  not  inherent 
in  the  model,  despite  its  simplistic  treatment  or  propaga¬ 
tion  and  other  possible  geometrical  effects. 

Reference  9  treats  the  case  of  nuclear  transitions,  for 
which  the  wavelength  is  smaller  than  the  atomic  spacing, 
but  '  >e  paper's  discussion  of  superfluorescence  is  less  in¬ 
formative  than  that  of  BL.  This  is  partly  because  the  BL 
model  uses  a  representation  for  the  collective  excitation 
operator  that  relates  it  to  the  lattice  structure,  whereas 
Ref.  9  uses  a  representation  that  relates  the  operator  to 
the  phases  of  the  emitted  photons  at  the  lattice  sites, 
thereby  missing  the  effect  of  the  regular  lattice  structure 
on  the  subsequent  radiation. 

Reference  10  gives  an  analysis  of  superradiar.ee  that, 
like  the  BL  model,  relates  the  collective  excitation  to  the 
lattice  structure  and  should  therefore  be  capable  of  yield¬ 
ing  the  short-wavelength  superfluorescence  conditions. 
Because  its  treatment  includes  more  geometrical  detail 
than  the  BL  treatment  for  the  long-wavelength  case  Ref. 
10  is,  in  fact,  able  to  relate  the  spatial  pattern  of  the  radi¬ 
ation  emitted  by  the  collective  atomic  array  to  the 
characteristic  shape  of  the  active  volume.  The  BL 
analysis  however,  is  restricted  to  the  case  in  which  the  ac¬ 
tive  volume  is  a  thin  needle. 

On  the  other  hand.  Ref.  10  shows  that  the  needle  shape 
provides  the  maximum  gain,  albeit  in  a  single  direction. 
Thus,  to  the  extent  that  it  is  possible  to  implement,  the 
needle  is  presumably  the  preferred  shape  for  a  superradi¬ 
ant  structure.  Therefore,  it  seems  useful  to  consider  the 
consequences  of  the  BL  model  for  the  short-wavelength 
radiation  that  characterizes  the  nuclear  case. 

The  interaction  Hamiltonian  originally  introduced  by 
BL  has  the  form 

H!lt)m~£v  2  2*k  0 k exp[/ ( UQ-oij  If  - ik-x;  ] 
v  y  k 


*/</>«  -XL  2  2*k(«^  “(a)/*(k  — o.i)  — H.c.  ]  . 

VV  •  k 

(4) 

where 

1  ^ 

/<*/>« TT  2  e*p(*iJ-*y>e*p[,‘(«y-‘»o> 0  •  (5) 

N  ;“i 

The  time-dependent  factors  in  (5)  are  separated  out  on 
the  assumption  that  they  can  be  replaced  by  ah  average 
over  all  atoms  and  that,  because  the  frequencies  e>;  asso¬ 
ciated  with  the  individual  atoms  are  uncorrelated,  they 
can  be  removed  from  /(ij,f)  and  subsumed  in  the  cou¬ 
pling  constant  gk,  which  then  becomes  time  dependent. 
This  time  dependence  of  gk  is  a  way  of  introducing  inho¬ 
mogeneous  broadening. 

At  this  point,  BL  introduce  the  assumption  that  the 
atomic  spacing,  given  by 

d=L/N  , 

is  much  smaller  than  the  photon  wavelength,  which  im¬ 
plies  that  the  function  /(q)  (wherein  the  time  depen¬ 
dence  is  no  longer  indicated  for  the  reason  just  noted)  is 
large  when  ij  vanishes  and  is  otherwise  small.  However, 
if  this  assumption  is  not  made  then  /<ij)  will  be  large 
whenever  ij-xJ*2trv;,  v;“0,±l,...  and  otherwise 
small. 

The  interaction  that  results  in  photon  emission  when 
cooperative  deexcitation  of  the  atoms  takes  place  can 
only  occur  when  /(k— a)  in  (4)  is  not  negligible  for  some 
value  of  k  and  a,  i.e.,  when 

(k— cr)*xy“2rvj  ,  (6a) 

where  vy  is  in  the  set:  0±1 . which  implies  that 

d /X-n/A'-v,  n=0,l . AT-1  ,  (6b) 


“H.c.  J  ,  (1) 

where  k  is  the  emitted  photon  wave  vector  for  a  particu¬ 
lar  mode,  Xj  is  the  position  vector  of  the  Jth  lattice  site, 
a  k  is  the  corresponding  photon  creation  operator,  is  a 
mean  reference  frequency,  and  r~  is  the  deexcitation 
operator  for  the  atom  at  the  Jth  lattice  site.  The  corre¬ 
sponding  excitation  operation  is  r*.  BL  also  define  the 
reciprocal-lattice  vectors  a,  each  component  of  which 
has  the  form  (2ir/Lln,  n  ■•0,1, . . . , N  —  1,  in  which  L  is 
the  corresponding  lattice  dimension.  They  define  the  col¬ 
lective  excitation  and  deex citation  operators  R±(a)  as 
functions  of  a  by 

it 

K  *(«)■»  2  rj±exP( ±ia-Xj ) ,  (2) 

/-i 

from  which  it  can  be  shown  that 

*f"'jr2*±(a)exP(T,'a*V  •  (3) 


where  v  is  in  the  set:  0,±1 .  In  (6)  it  is  tacitly  as¬ 

sumed  that  the  propagation  direction  is  along  a  narrow, 
needle-shaped  active  volume,  in  which  only  a  single¬ 
plane-wave  photon  mode  is  supported.  More  generally,  X 
should  be  interpreted  as  X^,  the  photon  wavelength  divid¬ 
ed  by  the  component  of  the  unit  vector  in  the  propaga¬ 
tion  direction  along  the  lattice  direction  defined  by  the 


It  follows  from  6(b)  that  the  possible  values  of  v  must 


be  non-negative.  Thus, 


kj-LANvj+nj),  0<nj<N,  vj>0  .  (7) 

On  setting  v;*0,l, ...  in  succession  it  follows  from 
(7)  that  for  Vj=Q,kj>d,  for  Vj*l,d >Xjid/2,  for 
vj-2,d/2>kj>d /3,  etc.  That  is,  for  any  given  photon 
wavelength,  one  and  only  one  value  of  v  exists  that  will 
satisfy  the  required  condition. 

Reference  9  gives  as  the  condition  for  a  Borrmann 
mode  the  relation 


lk-«|-lk| 


(8) 


On  substituting  (3)  into  (1)  it  follows  that 


for  some  reciprocal-lattice  vector  a.  To  satisfy  (8)  it  is 
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sufficient  that  k**e'  for  tome  reciprocal -lattice  vector11 
o’,  a  condition  that  (7)  implies.  That  is,  the  condition  im¬ 
plied  by  the  BL  model  for  supemdiance  in  a  crystal  is 
sufficient  for  coupling  to  a  Bomnann  mode,  the  propaga¬ 
tion  direction  of  which  will  be  determined  by  the  ratio  of 
the  wavelength  X  to  the  atomic  spacing  d  in  accordance 
with  (7).  . 


m.  VALIDITY  of  the  SEMICLASSICAL 
APPROXIMATION 

BL  show  that  their  model  implies  two  basic  conserva¬ 
tion  laws.  One  preserves  the  balance  between  emitted  ra¬ 
diation  and  stored  energy,  and  the  other  preserves  the 
Dicke  cooperation  eigenvalue  defined  in  terms  of  the 
time-varying  atomic  dipole  polarization  and  population 
inversion  states. 

The  first  law  has  the  form 

•7-  2  {<  («)>(/)]  + (Bj  >(/) 

dt  •“*»-** 

•-2K  2  [<><f(a)A(a))(r)]  ,  (9) 

— Sr"*. 

where  A  <  kg)  is  the  resonant  mode  of  the  internal  field,  ko 
is  the  vector  wave  number  at  resonance.  R  3  is  the  popula¬ 
tion  inversion,  and  K  (given  by  c/2 L,  where  L  is  the  axial 
length  of  the  active  volume  and  c  is  the  velocity  of  light) 
is  the  reciprocal  maximum  round  trip  transit  time  of  pho¬ 
tons  in  the  active  volume.  The  brackets  (  )  refer,  as  usu¬ 


where  g0  is  the  coupling  constant  at  resonance,  0  is  the 
volume  of  the  active  region,  and  TJ  is  a  time  constant 
due  primarily  to  inhomogeneous  line  broadening.  The 
unknown  quantities  appearing  in  Eq.  (11)  and  the  two 
conservation  laws  (9)  and  (10)  are  the  photon  number  ex¬ 
pectation  ( A  A),  the  atomic  inversion  expectation 
(R^),  the  photon  number/atomic  inversion  correlation 
(A  ARs),mm}  the  fluctuations  (R  *R  (R  }>  of  the 
atomic  dipole/inversion  vector  components. 

Equations  (9),  (10),  and  (11)  derived  from  the  BL  model 
do  not  form  a  complete  set  of  relations  for  all  of  the  ex¬ 
plicitly  involved  quantities  that  must  be  taken  as  indepen¬ 
dent  in  a  quantum-mechanical  treatment.  However,  in 
the  classical  approximation  the  relations 

<*!>-<*,>* , 

hold,  and  they  reduce  the  number  of  unknowns  to  3. 
Then  (9),  (10),  and  (II)  reduce  to  a  differential  equation, 
similar  to  that  derived  from  classical  mechanics  for  the 
•notion  of  a  pendulum,  and  a  corresponding  energy  rela¬ 


al,  to  the  expectation  value  of  the  operator  that  they  en¬ 
close,  and  parentheses  (  )  to  a  functional  dependence  on 
the  independent  variable  that  they  enclose.  The  second 
law,  in  which  R  *  and  R  ~  are  collective  dipole  moment 
operators,  has  the  form 

•7-  2  [(B  *(a)R  ~(a)>(r)] 

+  <Rj >(r)— <Aj)(f)  -0  .  (10) 

Any  Hamiltonian  system  would  imply  the  first  law 
The  second  is  analogous  to  and  formally  identical  with 
the  standard  conservation  of  angular  momentum  (result¬ 
ing  in  this  case  from  a  collection  of  pure  spin  states)  when 
the  total  angular  momentum  is  identified  with  the  Dicke 
cooperation  eigenvalue  associated  with  a  collective  tout 
angular  momentum  (spin)  operator  R,  and  the  angular 
momentum  vector  components  are  identified  with  Bj  and 
the  real  and  imaginary  parts  of  R  *. 

The  conservation  laws  are  therefore  physically  reason¬ 
able  in  their  own  right.  In  fact,  they  appear  to  be  quite 
general  and  could  be  regarded  as  essential  requirements 
for  any  model  based  on  the  collective  behavior  of  identi¬ 
cal  two-state  atoms. 

From  their  general  master  equation,  specialized  to  the 
case  of  two  identical,  independent,  single-resonant  modes. 
BL  also  derive  another  equation  involving  the  expecta¬ 
tion  values  of  the  atomic  and  electromagnetic  field  opera¬ 
tors: 


[<B*(a)B  ~(a>)(r>  +  2<  /4,<a)^(a)J{j)lt)]  ,  (Hi 


tion  which,  together,  do  form  a  dosed  system: 

2  <  A,(a)d(a))(t)=— ^j[^(r)]2e,/ri  , 

•  -ko.-k,  *g  0 

.  (12) 

1  Bo N  -t/T* 

d(f)+  K  +  — —  Ht) - e  Jsin^(f)=0, 

r;  » 

where  ^(  t)  is  a  modified  Bloch  angle  defined  by 

<B3X/)*[l+tfcos*(»)]/2 .  (13) 

With  relations  (12)  and  (13),  to  calculate  the  population 
inversion  and  the  emitted  electromagnetic  radiation, 
given  by 

<141 

2*o 

as  functions  of  time,  is  comparatively  simple. 

Thus,  where  the  classical  decorrelations  are  valid,  the 
BL  model  gives  a  firm  footing  to  a  numerically  tractable 
differentia)  equation,  permitting  straightforward  calcula¬ 
tions  of  the  most  important  quantities  associated  with 
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superfluorescence.  The  equation  includes  line  broadening 
and  accounts  fully  for  the  effects  of  atomic  stimulation  by 
the  local  field,  which  is  due  to  spontaneous  emission  by 
the  initially  excited  atoms. 

Further  on,  this  section  will  present  some  numerical  re¬ 
sults  obtained  from  the  BL  model  with  (he  aid  of  approx¬ 
imations  from4'1  tha,  unlike  the  classical,  preserve  first- 
order  quantum-mechanical  effects.  Those  results  indicate 
the  presence  of  large  quantum  fluctuations  during  the 
time  period  when  most  of  the  radiant  pulse  energy  is 
emitted.  This  is  somewhat  disturbing  because  the  validi¬ 
ty  of  the  classical  approximation  over  any  time  interval 
should  depend  on  quantum  fluctuations  being  small 
enough  to  be  neglected  during  the  interval. 

In  Ref.  12  Bonifacio,  Schewendiman,  and  Haake  report 
a  similar  finding  derived  from  an  earlier,  more  primitive 
version  of  the  BL  model.13  one  that  does  not  include 
stimulation  effects.  Physically,  the  earlier  model  (which 
is  Markovian  and  is  a  limiting  form  of  the  more  sophisti¬ 
cated  version  as  the  ratio  of  the  cooperation  time  to  the 
photon  propagation  time  becomes  large)  differs  from  that 
of  Ref.  5  by  virtue  of  the  fact  that  the  radiated  photons 
leave  the  active  volume  before  they  can  interact  with  the 
atoms.  As  a  result,  they  follow  the  atomic  state  changes 
adiabatically  and  do  not  produce  a  ringing  effect  in  the 
emitted  pulse. 

With  the  earlier  model,  the  authors  are  able  to  show  by 
direct  calculation13  that  quantum  fluctuations  are  large 
when  the  atomic  system  is  totally  inverted  initially.  But 
when  the  active  population  is  sufficiently  large  and  initial¬ 
ly  less  than  totally  inverted,  their  calculations  show  that 
the  fluctuations  are  small.  They  also  demonstrate  that, 
consistent  with  this  result,  the  classical  approach  is  valid 
whenever  the  initial  atomic  system  is  less  than  totally  in¬ 
verted. 

Although  the  results  in  Ref.  12  appear  to  validate  the 
use  of  the  classical  approximation  to  predict  spontaneous 
cooperative  radiation  whenever  the  atoms  are  not  totally 
inverted  initially,  the  more  sophisticated  Ref.  5  model 
does  not  necessarily  lead  to  the  same  conclusion  when 
non-Markovian  effects  are  important.  This  puts  the  va¬ 
lidity  of  the  semiclassical  approach  in  doubt  during  the 

_ I 


time  period  when  most  of  the  radiation  takes  place. 

Unfortunately,  the  region  of  validity  of  the  Ref.  5  ap¬ 
proximations,  which  take  into  account  quantum- 
mechanical  effects,  is  itself  uncertain,  at  least  for  the  case 
in  which  the  process  is  non-Markovian.  Thus,  calcula¬ 
tions  based  on  those  approximations  cannot  be  used 
directly  to  assess  the  accuracy  of  the  semiclassical  ap¬ 
proach  over  the  questionable  time  period. 

To  take  into  account  quantum  fluctuations  in  the  non- 
Markovian  case,  BL  make  two  approximations.  One  is 
the  Bom  approximation  which,  without  some  additional 
step  such  invoking  the  semiclassical  decorrelations,  does 
not  lead  directly  to  a  closed  system  of  equations  for  ex¬ 
pectation  values  of  photon  and  atomic  operators. 

BL  introduce  the  Dicke  states  [n,m  ),  which  are  eigen¬ 
vectors  of  the  photon  number  operator  A 1 A  and  the  in¬ 
version  operator  R  Jf  so  that 

A* A  |n,m  )“n|n,m  )  , 

>  “m \n,m  >  . 

In  terms  of  an  operator  ♦(»),  which  is  a  projection  of  the 
density  operator  of  the  system  that  leaves  the  atomic  ob¬ 
servables  unaltered,  BL  define  the  occupation  probabili¬ 
ties 

p(m,f)*=2  (n,m  id>U)|n,m  >  , 

■ 

the  photon  number  expectations 
Anm.f)*^  <n,m|d>(/)|»i,m  )  , 

and  the  quantities 

Llm.O-ReJfln  +lMn  +2)],/2 

ft 

X  <n +2,m  — 2|d>(r)|n,m  )  . 

In  the  Dicke  sute  representation,  the  Bora  approxima¬ 
tion  leads  at  first  to  a  finite  system  of  integrodifferential 
equations  for  the  occupation  probabilities  p(m,t)  of  the 
Dicke  sute  basis  vectors  |r,m  ): 


.  2go  ri.  -h+si/itT ~Kn~t> 

p{m,t)= - I  ds  e  1 

v  Jo 

X\g[m)p{m,s)—g{m  +1  )p(m  +l,s)+[g{m)+glm  +  l))W(m,s) 

—g(m  +lW(m  +l,i)-|lmW(i«  —  l,s)+g  ,/l(m)g  ,/3(m  —  UL(m,j) 

—2g,/3(m)g,/2(m  +l)Hm  +  l,s)+gl/7lm  +l)gl/1lm  +  2)Z.(m  +2,s)j  ,  (15) 

where 

<jAr’+m)<jAr-m  +1)  for  -  $N' ’  <  m  <  , 

0,  otherwise, 

and  N'  is  the  number  of  emitters.  Unfortunately,  the  number  of  unknown  quantities  to  be  determined  from  Eq.  (15)  is 
larger  than  the  number  of  equations  in  the  system.  When  the  authors  enlarge  the  system  to  include  time  derivatives  of 
the  unknowns  other  than  the  Dicke  state  occupation  probabilities,  the  resulting  set  of  equations  is  still  not  closed  be¬ 
cause  it  contains  new  unknowns,  consisting  of  higher-order  moments  of  the  photon  and  atomic  operators.  Obviously, 
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repeating  the  process  indefinitely  will  result  in  an  infinite  hierarchy  of  (finite)  systems  of  equations,  involving  moments 
of  ever  increasing  order. 

BL  observe  that  a  simple  way  of  getting  a  closed  system  is  to  drop  the  photon  number  expectation  values  (but  retain 
their  derivatives  on  the  left-hand  side)  and  all  second-order  moments  from  the  first  set  of  equations  in  this  hierarchy 
that  involve  time  derivatives  of  just  the  Dicke  state  occupation  probabilities  and  the  photon  number  expectations. 
However,  that  procedure  should  be  valid  only  under  conditions  that  would  justify  substituting  the  earlier  Markovian 
model.  It  would  therefore  add  nothing  new,  serving  only  to  verify  that  their  earlier  model  is  a  limiting  case  of  the  more 
general  non-Markovian  model 

Their  next  step,  then,  is  to  keep  the  equations  foi  which  the  left-hand  side  involves  time  derivatives  of  the  Dicke  sute 
occupation  probabilities  and  the  photon  expectation  values  and  retain  the  Dicke  state  occupation  probabilities  and  the 
photon  expectation  values  on  the  right-hand  side,  but  to  drop  all  second-order  and  higher  moments.  The  result  is  a 
larger,  but  still  closed,  system: 


e 


X  |g(m)p(m,s)— g(m  +  l)p(m  +  l,s)  +  [g(m)+g(m  +  l)jAf(m,s) 
— g(m  +1  )N(m  +1  ,s)—g(m)AT(m  — l,j))  , 


Mm.O—  —  2JfAf(m,0+ -  f  ds  t  1 

v  Jo 


(16a) 


X  |g(m  +l)[p<«  +  l,r)+Af(i»i  +  l,s)— Af  (m,s)]j  . 


(16b) 


Significantly,  the  equations  resulting  from  this  approxi¬ 
mation  imply  a  relationship  between  atomic  and  photon 
operator  expectation  values  that  is  also  implied  by  the  ex¬ 
act  (i.e.,  with  no  approximations)  operator  equations  of 
the  model.  Therefore,  the  three  previously  discussed 
equations,  (9),  (10),  and  (11),  follow  from  (16).  In  addi¬ 
tion.  (16)  guarantees  conservation  of  probability:  The 
sum  of  the  Dicke  sute  occupation  probabilities  satisfying 
(16)  must  remain  constant  over  time.  This  result  follows 
from  the  idenuty 

2  f(m  +  l)F(m  +  1)—  2  glm)F(m) , 

in  which  F(m)  can  be  any  function. 

Quantum-mechanical  fluctuations  are  responsible  for 
the  triggering  of  the  superfluorescent  pulse  from  the  in¬ 
verted  population.  It  has  been  argued  that  when  these 
fluctuations  are  small,  classical  solutions  give  good  ap¬ 
proximations  to  the  superfluorescent  pulse  emission  and 
that  these  fluctuations  are  small  at  the  maximum  emis¬ 
sion  rate.u_l5 

To  check  this  assumption  and  study  the  characteristics 
of  solutions  obuined  from  Eq.  (16)  we  calculated  the  em¬ 
itter  pulse  given  by 

7(0-2  (H> 

m 

and  the  variance  defined  as 

=  .  (18) 
_  \  m  I 


FIG.  1.  SF  pulse  intensity  (a)  and  variance  (b)  calculated  for 
the  system  with  N  — 100  (  —  SOS m  £50)  and  AT  — 10.  The  ini¬ 
tial  condition  is  the  totally  inverted  state  with  the  probability 
Pim,0)m0  for  all  m  except  for  m  “  50  for  which  F(50,0)“  1. 
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For  the  calcu'itions  we  will  discuss  here  we  assumed 
100  panicles  so  that  N  «*  100  and  m  varies  from  —SO  to 
+50.  We  assumed  different  initial  conditions  for 
the  probabilities  P(m,t)  and  values  of  K  to  obtain 
both  pure  and  oscillatory  SF.14  In  our  calculations  we 
measured  time  in  units  of  the  cooperation  time 
r,  “( 1  /g0  K  V/N)xn.  Figure  1  shows  the  result  of  our 
calculations  undo1  pure  SF  conditions  with  K  “  10.  The 
initial  probabilities  are  f(n,r)"0  for  all  m  except 
m  “50,  which  give  a  fully  invened  state.  Note  that  the 
variance  is  maximum  when  the  pulse  is  being  emitted 
[compare  Figs.  1(a)  and  1(b)].  These  results  indicate  that 
the  quantum  fluctuations  are  not  necessarily  negligible 
during  the  course  of  the  emission  process.  Also  note  in 
Fig.  2  the  shift  of  the  state  population  from  the  inverted 
state  [right  side  in  (a)]  at  early  times  to  the  broader  distri¬ 
bution  at  intermediate  times  (d),  (e),  and  (f)  and  finally  the 
ground-state  population  at  higher  times  (g)  and  (h). 

When  we  attempt  to  run  the  same  calculation  with 
K  <2^2  an  instability  develops:  The  intensity  oscillates 
wildly  until  the  calculation  blows  up.  The  violent  osdlla- 


(•) 


FIG.  2.  Probability  distribution,  P(m,i)  at  different  times  for 
the  states  —505m  <30,  used  in  the  calculation  of  the  SF  inten¬ 
sity  and  variance  shown  in  Fig.  1.  The  times  at  which  the  prob¬ 
ability  Plm.t 1  is  shown  are  I  “0.05,  4,  g,  10,  13,  IS,  and  24  in 
(a)-(h),  respectively. 


(b) 


FIG.  3.  Radiative  poise  and  quantum  fluctuations  (variance)  for  different  K  values  which  are  0.0269  units  in  (a),  0.269  in  (b),  0.8  in 
(cl,  and  2.69  in  (d).  The  initial  distribution  for  (ml  is  Pouaon  with  P,  ■  10~’  peaked  at  m  —  - 50 and  N(m)“0. 
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tions  and  blow  up  can  be  reduced  by  starting  with  the 
population  closer  to  the  ground  state,  as  shown  in  Fig.  3, 
where  the  Poisson  distribution  is  peaked  at  the  ground 
state  M  “  —  50.  In  Fig.  3  for  K  values  from  0.0269  to 
2.69,  below  the  critical  value  of  2^2  the  emitted  pulses 
remain  stable  and  the  variance  is  positive.  This  is  not  the 
shown  in  Fig.  4  when  the  initil  distribution  is 
peaked  at  m  *—46  [P(m,r)*l  for  m  *—46  and  zero 
otherwise]  and  K  *0.0269,  the  intensity  goes  slightly 
negative  and  the  variance  oscillates  between  positive  and 
negative  values.  With  the  same  initial  distribution  but 
with  higher  K  *2.69  (still  less  than  2v/2)  the  tingle  pulse 


appears  and  the  variance  stays  positive  as  shown  in  Fig. 
3.  In  Sec.  IV  we  examine  the  source  of  the  instabilities 
and  the  conditions  for  their  occurrence.  We  also  discuss 
a  closed-form  analytical  solution  obtained  for  a  system 
with  two  particles. 


IV.  NUMERICAL  INSTABILITY 
IN  THE  BL  APPROXIMATION 

Because  the  approximate  equations  (16)  yield  the  exact 
basic  conservation  laws  of  the  BL  model,  it  might  be  ex¬ 


pected  that  their  solutions  would  be  physically  well 
behaved.  However,  numerical  calculations  indicate  oth¬ 
erwise. 

After  a  certain  time  interval,  before  the  emitted  pulse 
reaches  its  maximum  amplitude,  quantities  derived  from 
the  approximate  solution  become  nonphysical  in  at  least 
two  respects.  First,  although  conservation  of  total  proba¬ 
bility  is  still  satisfied,  individual  Dicke  state  occupation 
probabilities  become  negative.  Second,  although  conser¬ 
vation  of  energy  is  still  satisfied,  individual  Dicke  state 
photon  number  expectation  values  also  become  negative. 

The  case  in  which  the  number  of  atoms  in  the  active 
volume  is  limited  to  two  is  simple  enough  to  be  treated  in 
detail  analytically  if  inhomogeneous  line  broadening  is 
neglected.  An  investigation  of  it  using  the  Laplace  trans¬ 
form  reveals  that  all  states  below  the  maximum  Dicke  oc¬ 
cupation  number  (which  is  two  in  the  case  considered)  ex¬ 
hibit  a  resonance  phenomenon;  i.e.,  some  of  the 
transformed  solution  functions  have  double  poles.  In  the 
inverse  transform  domain  such  a  function,  which  would 
otherwise  be  exponential  or  trigonometric,  must  have  the 
time  variable  as  a  factor. 

In  the  absence  of  double  poles,  when  parameters  are 


FIG.  3.  ( Continued ). 
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0  12.0 
time  (units  of  te) 

FIG.  4.  Pulse  (a)  and  variance  (b)  for  (be  initial  conditions 
p(3,0)-l,p(m.0)-0.  m*5,  A-lOO,  and  K-0.0269. 


time  (units  of  te) 


FIG.  S.  Poise  (a)  and  variance  (b)  for  the  initial  conditions 
*(5,0)-  l,*(m,0>-0,  N  « 100,  and  K  -2.69. 


such  that  ringing  does  not  occur,  the  solution  functions 
decay  exponentially  with  time.  When  the  converse  is 
true,  i.e.,  parameters  are  such  that  ringing  does  occur, 
the  solution  functions  have  factors  that  are  linear  com* 
binations  of  trigonometric  functions  of  time.  At  least  in 
the  second  case,  the  additional  time  variable  factor  im¬ 
posed  by  a  double  pole  guarantees  increasingly  larger  os¬ 
cillations  with  linearly  increasing  amplitudes  that  must 
eventually  produce  negative  expectation  values. 

In  principle,  a  similar  analysis  for  an  arbitrary  number 
of  atoms  could  be  carried  out  in  the  same  way,  using  the 
Laplace  transform,  since  the  equations  can  be  treated  re¬ 
cursively  in  pairs,  no  matter  what  the  tout  number  of 
atoms  may  be.  Because  of  the  symmetry  of  the  inversion 
operator  eigenvalues  about  zero,  the  coefficients  due  to 
the  lower  Dicke  states  will  always  produce  double  poles 
for  the  transformed  solution  functions,  just  as  in  the  case 
of  two  atoms. 

In  fact,  in  Ref.  12  the  authors  make  the  same  observa¬ 
tion  in  connection  with  the  earlier  Markovian  model,  i.e., 
that  due  to  this  coefficient  symmetry,  double  poles  must 
always  occur  in  the  transformed  solution  functions. 
However,  in  Ref.  12  the  remark  is  made  in  passing, 
without  noting  the  consequence  that,  in  approximate 
solutions  for  a  more  general  model,  such  resonances  may 
lead  to  physically  impossible  negative  expectation  values. 

Equations  (16)  arc  approximate  relations  for  the  occu¬ 
pation  probabilities  p(m.t)  of  the  atomic  Dicke  states 
|r,m  land  the  corresponding  photon  expectation  values 
iY(m.t)  of  the  radiating  electromagnetic  field,  given  a 
particular  value  iN'/2){N'/2+ 1)  for  the  atomic 
cooperation  eigenvalue.  For  the  sake  of  an  analytic  ex¬ 
ample  we  now  take  the  “inhomogeneous  broadening”  de¬ 
cay  time  r?  to  be  infinite.  Time  is  measured  in  units  of 
the  cooperation  time  re;  thus,  the  quantity  2g§/i>  be¬ 
comes  4  /N'. 

Along  with  Eq.  (16)  the  standard  initial  conditions 


N(m, 0)-0  , 


*(m,0)-0, 


(19) 


will  be  assumed  here.  These  conditions  imply  that  initial¬ 
ly  all  N'  atoms  are  excited  and  the  photon  field  is  in  the 
vacuum  state.  Together  with  (16)  they  imply  that 
N(N'/2,t)=0. 

Also,  with  neglect  of  TJ,  the  solution  of  (16)  can  be 
found  by  means  of  the  Laplace  transform.  Then,  setting 
pimfz)  and  N{m,z)  for  the  Laplace  transforms  of p(m,t) 
and  (16)  is  equivalent  to 


r+2*  + 


4 

z+K 


*0  , 


z+K 


N 


4» 

,+  * 

P 

N' 

2  ' 

^  z+K 

T’2 

for  the  highest  m  values  and 
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*+2JC  +  ~p  Af(m  — 1.2} - ~—p(m,x) 

X  TA  X  T A 


r N(m,z)  , 


~^—S(m  - 1,2)+  *  +  ~7r 

X  "»  A  x  TA 


r,+r,*i 

z+* 


AT(m.z) 


which  in  terau  of  the  present  notation  would  be 

Jf-1. 

IT  oscillation  at  a  natural  frequency  u  occurs  for  some 
value  m'  of  m,  then,  since  the  corresponding 
and  p(m'.t)  will  contribute  to  the  sources  in  the  equa¬ 
tions  for  Nlm'-2,t)  and />(«'—  l,f),  each  of  those  quan¬ 
tities  will  also  have  a  term  oscillating  with  the 
frequency  a>,  etc.  Thus,  if  there  is  a  natural  frequency 
again  equal  to  a>  for  some  m  on',  the  corresponding 
S{m  —  1  ,t),p(m,t)  will  exhibit  resonant  vibrations. 


Ym  + 1  -  Ym  + 1  ————— 

- —-—Sim  +  l,z) - —--pirn  +  l.z) 

r+JT  z+K 

for  the  rest,  where 

,  Alai .  an 

Ym  S' 

The  quantities  Sim  —  l,z)  are  naturally  matched  with 
the  quantities  p(m,z)  in  Eq.  (20)  in  the  recursive  pairs, 
the  quantity  SiS'/2,t)  already  having  been  found  to  be 
zero  identically.  The  determinant,  Det(z),  of  the  matrix 
on  the  left  side  of  Eq.  (20)  is  given  by 

_  .  .  _ _ f  m  Ym  Ym * 

Det(r)=  z+2K  — z  +  r 

z  +K  x+KT  iz+Kr 
**z2+2Kz  +2 Ym  .  (22) 


0  time  (unite  of  9 


£  .1 


Dm#  (unite  of 


solve  Eq.  (20).  is  given  by 

— 1 W  it  0.7 

91^  A 

zJ+ifz+r- 

Ym 

p  f\ 

Det(z) 

Det(z) 

r  \ 

A  ~,*B 

Ym 

z3+3JCz+21C1r. 

(23) 

t  \ 

Det(z)  * 

Dci(z) 

S  o  f  V _ 

The  zeros  of  Det(z)  determine  whether  the  solution 
pair  Nim.t)  and  p(m,f)  for  a  given  value  of  m  oscillate 
with  a  natural  frequency  that  does  not  depend  on  the 
right  side  of  Eq.  (20).  If  the  zeros  are  both  real,  those 
quantities  do  not  have  a  natural  oscillation  frequency;  if 
the  zeros  are  complex  (in  which  case  they  occur  in  conju¬ 
gate  pairs)  they  have  a  natural  oscillation  frequency  o> 
given,  except  for  sign  which  is  conventionally  positive,  by 
the  imaginary  part  of  either  zero. 

The  zeros  are  given  by 

z - -K±VK1-lYm  -  -K±V~K*-ig{m)/S' .  (24) 

The  condition  for  a  natural  oscillation  frequency  corre¬ 
sponding  to  a  given  m  is  therefore 

K<2Vlg(m)/S'  .  (25) 

The  definition  of  the  g(m),  along  with  (25),  together 
imply  that  ringing  occurs  in  the  system  if  and  only  if*7 

K<  2^2  . 

This  is  a  more  precise  condition  for  the  existence  of  a 
non-Markovian  stimulation  effect  than  the  Ref.  5  condi- 
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time  (unite  of 


0  time  (units  of  x^  9 


0  time  (unite  of  tfi)  9 

FIG.  6.  SF  pulses  given  by  /(f)* ^^Sim.i)  for  S'm 2  cal¬ 
culated  from  Eq.  (32).  In  (a)  K  -0.03,  (b)  K  -0.3.  (c)  A  -3.  (d 
K -3.  (e)  K  -JO. 
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It  is  evident  from  (24)  that  this  can  occur  if  and  only  if 
g{m‘)mg(m)  . 


g  — — ft  “(m+1  KN'-p)  , 


if  and  only  if 
tf’*M+v+ 1 


N'-2v+l,  «8) 

so  that  to  avoid  resonance  and  still  have  oscillation,  N' 
must  be  odd. 

Also,  on  comparing  Eq.  (26)  with  Eq.  (28),  it  is  found 
that  this  can  only  occur  when  m  for  which  value 


g(m)«g  - 


Then,  the  natural  frequency  oscillation  condition  Eq.  (25) 
becomes 


*<<Ar+l)v/27N5~v'23r  . 


.  .  To  avoid  resonance  due  to  the  natural  frequency  for  a 

It  follows  that  for  any  JV>1  and  any  v  such  that  j,  value  of  m,  the  condition 

0<  v<JV'  — l,  *  ___ 

Jlf>2v/g(j)/^'“2v'(jV'  +  3KN;-l)/2N'  (30) 

t  v+1-^f  •  (26)  must  be  satisfied. 

In  summary,  the  only  case  in  which  there  e  oscil- 

It  follows  from  Eq.  (26)  that,  if  a  natural  frequency  ^0“  the  *y»«m  but  no  resonance  occurs  cn 

occurs  for  some  value  of  m  given  by  S‘/2-v  or  v'2(JV'-+-3KJV,'~ I )>JV*<JC  kW+DvTTs'  .  (31) 

v+ 1  —S'/ 2,  resonance  must  occur  for  some  lower  value 

unless  The  case  of  S'*=  2  provides  an  example  for  which  the 

v,  v,  explicit  analytical  resulu  are  reasonably  simple.  On  ap- 


T_v.v+1__. 

Condition  (27)  is  equivalent  to 


(27)  plying  the  matrix  inverse  Eq.  (23)  to  Eq.  (20)  and  invert* 
ing  the  Laplace  transforms  obtained  thereby,  the  follow¬ 
ing  solutions  are  obtained: 


A((l,r)»0,p(l,f)~— ^-je“*Tl— costVT— jpr  4-^)] , 

O  "“A 


fiVgTF  .  *J~4 
j - .  cos#*  -  , 

JV(-l,rl- — ^^*-*'[1— cosv/g— *ir-v'(8-/:J)/2rsinV''8— K^r] , 
p{0,t)~^Si-\,t)+  '/(8-/rr)/2sinV'8-/:Jf 


*=£■  t  'o.vtt*,*  (t-ffn,*vT=P,  , 


pi  —  l,r)«  1  —piO,t)—p (l,t) . 


The  last  equation  comes  from  the  conservation  of  proba¬ 
bility 

2p(m,r)-l  , 


show  some  calculated  resulu  from  Eq.  (32).  These  resulu 
compare  well  with  the  calculations  obtained  using  the 
previous  numerical  methods  of  solving  (16)  for  the  same 
input  parameters. 

V.  CONCLUSIONS 


which  follows  identically  from  Eq.  (16)  and  the  assumed  V.  CONCLUSIONS 

initial  conditions. 

The  emitted  pulse  is  given  by  the  sum  of  the  individual  The  BL  mean-field  model  of  superfluorescence  due  to 
photon  contributions  as  in  Eq.  (17).  Figures  6(ai-6(e)  an  array  of  identical  two-level  atoms  has  an  important 
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advantage  not  shared  by  other,  more  recent  and  more 
comprehensive  models.  It  deals  quantum  mechanically 
with  emitters  in  a  crystal  lattice,  which  is  the  most  likely 
structure  for  the  y-ray  laser. 

The  model  can  be  applied  to  a  similar  structure  com¬ 
posed  of  nuclear  emitters  that  radiate  electromagnetic  en¬ 
ergy  at  wavelengths  smaller  than  the  interatomic  spacing. 
It  predicts  that  superfluorescence  in  this  case  is  possible 
only  at  wavelengths  satisfying  conditions  that  are 
sufficient  for  Bornnann  coupling  to  occur.  However, 
since  the  BL  model  is  one  dimensional  it  cannot  distin¬ 
guish  among  the  various  Borrmann  modes  that  Trammel, 
Hutton,  and  Hannon"  have  discussed  in  connection  with 
a  three-dimensional  crystal. 

Solutions  of  the  BL  model  will  approach  a  classical 
limit  some  time  after  superfluorescent  radiation  begins; 
however,  in  some  cases  the  early  time  during  which  the 


classical  approximation  remains  invalid  lasts  until  after 
the  peak  of  the  radiated  pulse  has  emerged.  On  the  other 
hand,  it  is  found  that  in  some  cases  the  iteration  process 
used  to  approximate  the  quantum-mechanical  solution, 
which  is  required  during  the  early  stages  of  the  emission 
process  before  the  classical  limit  provides  a  valid  approxi¬ 
mation,  can,  itself,  be  invalid  because  of  a  numerical  in¬ 
stability  inherent  in  the  process. 
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Abstract 

A  theory  of  nuclear  superfluorescence  (SF)  based  on  the  Haake-Reiboid  model  for  the 
atomic  case  is  presented.  Certain  modifications  of  the  model  make  it  possible  to  take  into 
account  some  effects  that  are  more  important  in  nuclear  than  in  atomic  SF:  attenuation, 
competing  transitions,  finite  pumping  times,  and  both  homogeneous  and  inhomogeneous  line 
broadening.  Results  of  some  explicit  calculations  illustrate  the  influence  of  these  effects  on  the 
radiated  pulse.  Similar  calculations  are  used  to  examine  the  feasibility  of  observing  SF  using  the 
38.6  keV  transition  in  “Co. 


L  Introduction 

Superfluorescence  (SF)  (Refs.  1-4)  is  an  example  of  the  spontaneous,  cooperative  emis¬ 
sion  of  coherent  radiation  by  a  collection  of  identical  molecules,  atoms  or  nuclei  Intense, 
directed  pulses  of  duration  much  shorter  than  the  spontaneous  emission  (SE)  lifetime  of  an 
individual  radiator  characterize  the  radiation,  which  has  been  observed  experimentally  for  atoms 
and  molecules  (Ref.  2),  e.g„  in  CHjF,  HF,  Na,  Ca,  Tl  KC1,  Sr,  and  Li.  Existing  theoretical 
models  (Ref.  2)  have  successfully  explained  the  observations. 

However,  thus  far  no  one  has  reported  an  observation  of  nuclear  SF,  nor  have  the  theor¬ 
etical  models  used  to  analyze  the  corresponding  nuclear  phenomena  been  adequate.  This  is 
unfortunate  since  SF  (Refs.  3, 6)  will  probably  be  the  most  important  emission  process  in  a  y-ray 
laser  operation.  The  emission  time  of  a  cooperative  pulse  is  inversely  proportional  to  the 
number  of  cooperating  emitters,  and  is  therefore  likely  to  be  much  shorter  than  the  natural 
nuclear  lifetime.  Also,  the  peak  intensity  of  the  radiated  pulse  will  be  proportional  to  the  square 
of  the  number  of  cooperating  emitters  and  therefore  much  larger  than  the  intensity  of  the  inco¬ 
herent  SE  which  is  just  proportional  to  the  number  of  emitters. 

This  paper  takes  into  account  certain  experimentally  achievable  parameters  that  are  more 
important  for  nuclear  than  for  atomic  SF.  Its  theoretical  treatment  is  based  on  the  Maxwell- 
Bloch  equations  introduced  by  Haake  and  Reibold  (Ref.  7)  to  treat  superfluorescent  phenomena 
in  multi-level  systems,  but  modified  here  to  include  most  of  the  important  nuclear  effects.  We 
apply  the  theory  to  a  real  nucleus,  “Co,  created  by  thermal  neutron  irradiation  of  s*Co,  and  show 
under  what  conditions  nuclear  SF  is  possible  in  this  system. 
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II. 


Fig.  1  depicts  the  Haake-Reibold  model  for  SF,  including  some  parameters  that  cannot  be 
neglected  in  the  nuclear  case.  As  indicated  in  the  diagram  on  the  left,  level  4  decays  at  the  rate 
Y  in  a  pumping  process  that  populates  level  3.  The  inverted  population  of  level  3  then  decays 
with  the  rate  T  to  the  ground  state  at  level  1,  collectively  emitting  SF  radiation  that  is  proportion¬ 
al  to  the  square  of  the  population.  Appearing  on  the  right,  the  Maxwell-Bloch  equations  govern 
the  relation  between  the  populations  N«,  Nj,  N,  of  the  three  levels,  the  electric  field  E*,the  col¬ 
lective  polarization  Rs,  and  the  source  term  £*,  all  of  which  are  functions  of  distance  and  time.1 
In  general,  the  electric  field  satisfying  the  equations  is  a  pulse  propagating  from  left  to  right 


9N4 

•ar=-^4 

3N, 

—  (eV+e-r-)  -  r  N3+  tn4 

3N.  . 


=  +  (EV+E*R‘)  +  rN, 


^ = (N3-Nj)  et-  \  (T + r#)R±+  K 

^=g(t)RT-^iE±, 


The  unit  of  time  is 


8»t0 

TSF*3XV 

where  1  is  the  cavity  length,  X  is  the  wavelength  of  the  emitted  photon,  t0  is 
the  natural  radiative  lifetime,  and  p  is  the  density  of  cooperating  nuclei 


Figure  1 


The  Haake-Reibold  Model 


In  the  Maxwell-Bloch  equations  t  represents  retarded  time,  so  that  the  last  equation 
implies  an  electric  field  propagating  from  left  to  right  Also,  the  source  term  represents  noise 
due  to  fluctuations  of  the  vacuum  and  is  therefore  a  stochastic  quantity,  the  properties  of  which 
Ref.  8  derives  using  a  quantum  electrodynamic  argument 


’These  quantities  are  quantum  mechanical  operators,  but  because  the  source  term  becomes 
negligible  after  a  short  time  and  the  other  operators  all  commute  with  each  other,  they  can  be 
treated  as  c-numbers  and  identified  with  their  expected  values. 
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The  SF  model  depicted  in  Fig.  1  includes  the  following  phenomena,  all  of  which  are 
important  for  the  nuclear  case:  electromagnetic  attenuation,  natural  decay,  competing  transitions, 
homogeneous  and  inhomogeneous  line  broadening,  pumping  rate,  and  relaxation.  Various 
parameters  appearing  in  the  Maxwell-Bloch  equations  account  for  these  phenomena  and 
determine  their  ultimate  effect  on  the  system  radiation. 

As  indicated  by  the  first  equation,  the  reciprocal  of  y  is  the  effective  pumping  time.  The 
form  of  the  second  equation  shows  that  r  is  the  natural  decay  rate.  A  quantity  p  in  the  last 
equation  obviously  causes  attenuation  of  the  propagating  electric  field. 

The  homogeneous  line  broadening  mechanism  differs  from  that  of  the  inhomogeneous 
broadening.  The  fourth  equation  indicates  that  the  quantity  T*.  which  represents  dephasing 
associated  with  homogeneous  line  broadening,  simply  adds  to  the  effect  of  the  natural  decay. 
However,  the  function  g(t)  in  the  last  equation  is  responsible  for  taking  into  account  inhomoge¬ 
neous  line  broadening.  The  form  of  g(t)  depends  on  the  statistical  distribution  of  frequencies 
associated  with  random  dephasing  factors  contributed  at  a  particular  location  by  photons  emitted 
at  other  positions  Xj  along  the  active  region.  If  the  distribution  is  Lorentzian  g(t)  has  the  form 

qW  -  g0eTt- 

The  quantity  go  is  the  coupling  constant  relating  the  electric  field  to  the  polarization. 


m.  Predictions  of  the  Modified  Haake-Reibold  Model 

Using  a  Monte  Carlo  procedure  to  simulate  the  stochastic  source  numerical  solutions 
of  the  Maxwell-Bloch  equations  in  Fig.  1  lead  to  the  curves  given  in  Figs.  2-4,  which  illustrate  the 
effect  of  attenuation  and  line  broadening  on  the  radiated  SF  pulse.  Fig.  2  shows  how  the  peak 
intensity  and  the  delay  time  of  the  peak  vary  as  functions  of  p  for  a  cooperation  number  N  equal 
to  10s.  Fig.  3  illustrates  the  effect  of  homogeneous  and  inhomogeneous  line  broadening  on  the 
SF  pulse  shape,  using  curves  representing  intensity  as  a  function  of  time  for  this  purpose.  Fig.  4 
shows  the  delay  time  as  a  function  of  homogeneous  and  inhomogeneous  broadening.  The  figures 
use  a  parameter  a  equal  to  the  ratio  of  T,  to  the  natural  decay  rate  T  and  a  parameter  b  equal 
to  the  ratio  of  r*  to  r  as  convenient  measures  of  line  broadening. 

Ref.  9,  reporting  on  observations  of  combined  SF  and  amplified  spontaneous  emission 
(ASE)  produced  by  a  KC1:0'2  atomic  system,  gives  experimental  curves  showing  the  effect  of  line 
broadening  on  the  delay  time  of  emitted  pulses  of  both  types.  Calculations  based  on  the  Haake- 
Reibold  model  used  to  predict  SF  and  ASE  agree  quite  well  with  the  Ref.  9  results  when  reason¬ 
able  values  are  chosen  for  parameters  that  could  not  be  measured  or  otherwise  determined. 

Unfortunately,  no  one  has  yet  reported  a  similar  experiment  involving  nuclear  rather  than 
atomic  emissions.  However,  Ref.  10  gives  a  list  of  candidates  to  replace  KCliOj  in  such  an  ex¬ 
periment,  one  of  which  is  “Co  prepared  in  the  excited  isomeric  state  by  thermal  neutron  pum¬ 
ping  of  "Co. 

Aside  from  this  pumping  mechanism,  other  advantages  of  “Co  are:  (1)  its  lifetime,  906 
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sec,  which  would  allow  long  pumping  times  for  population  inversion  and  therefore  would  involve 
achievable  thermal  neutron  fluxes;  (2)  the  fact  that  its  natural  lifetime  is,  nevertheless,  relatively 
short;  (3)  its  relatively  good  Mdssbauer  effect;  (4)  its  low  losses  due  to  such  factors  as  internal 
conversion  and  electronic  attenuation. 


Figure  2 


Linear  Attenuation  Coefficient  p.  (cm*1) 


SF  Peak  Intensity  and  Delay  Time  as  Functions  of  Attenuation 


Figure  3 


SF  Pulse  Shape  Dependence  on  Line  Broadening 


Figure  4  SF  Delay  Time  Dependence  on  Line  Broadening 


Table  1  gives  the  characteristic  atomic  and  nuclear  parameters  associated  with  a  single 
“Co  nucleus.  It  omits  inhomogeneous  broadening  since  that  phenomenon  involves  multiple  nuc¬ 
lei  In  fact,  the  only  estimates  of  inhomogeneous  broadening  in  connection  with  a  long  lived, 
MOssbauer  effect  isotope  that  are  available  in  the  literature  are  for  107 Ag  (Ref.  11)  and  ,wAg 
(Refit.  12-14)  and  are  based  on  self-attenuation  measurements.  The  earliest  estimates  of  the  in¬ 
homogeneous  broadening  parameter  a  were  of  the  order  10s,  but  more  recent  ones  (Refs.  13, 14) 
are  much  smaller,  ranging  from  30  to  200. 

Table  1  gives  12  cm'1  for  the  electronic  attenuation  p.  However,  exploiting  the  Bomnann 
effect  in  a  single  crystal  can  reduce  that  value  to  -0.01  cm'1. 

In  a  generic  example.  Fig.  3  compares  the  “Co  SF  intensity  as  a  function  of  time,  calcul¬ 
ated  by  means  of  the  Haake-Reibold  model,  with  a  similar  calculation  of  the  SE.  It  shows  the 
number  of  photons  emitted  per  sec  per  nucleus  for  both  cases  separately  and  combined.  Whe¬ 
ther  the  SF  is  detectible  in  an  experiment  will  depend  on  whether  the  disturbance  it  causes  in  the 
combined  SF  and  SE  is  measurable. 

If  the  fluorescent  material  is  a  single,  acicular  shaped  crystal  the  SF  will  radiate  in  an 
endfire  mode  in  a  narrow  beam,  with  a  high  gain  relative  to  an  isotropic  radiation  pattern  charac¬ 
teristic  of  SE,  due  to  natural  decay.  If  the  diameter  d  of  the  active  region  is  related  to  its  length 
1  and  the  wavelength  1  by  the  Fresnel  condition 

fior  a  minimum  beamwidth  p,  Le.,  angular  divergence  due  to  diffraction,  the  beamwidth  is  given 
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Table  1 

Characteristic  Parameters  for  “Co 


Parameter 

Value 

Transition  energy 

58.6  keV 

Wavelength  on  resonance 

(UA(2  x  104  cm) 

Natural  lifetime 

628  sec 

RecoQless  fraction 

0304 

Linear  attenuation  coefficient 

12  cm*1 

Interval  conversion  coefficient 

483 

Particle  density  (in  solid  form) 

8.97  x  10s  cm4 

Thermal  neutron  Cross-section  J*Co+n  -  “Co 

20  barns 

Symbol 


E. 


Figure  5  Comparison  of  Spontaneous  and  SF  Photon  Emission  from  Inverted  “Co  System 


If  the  SF  beamwidth  is  5x10-*  rad,  which  corresponds  to  a  solid  angle  of  6.25  x  10-1®* 
sterad,  the  gain  relative  to  an  isotropic  polyciystaliine  radiator  will  be  -6.4  x  10*.  Neglecting 
inhomogeneous  broadening  and  attenuation.  Table  2  gives  the  normalized  photon  emission  rate 
for  both  the  case  of  a  polycrystalline  sample,  labelled  Experiment  A,  and  a  single  crystal,  labelled 
Experiment  B,  for  different  values  of  the  thermal  neutron  Qux  used  to  pump  the  **Co  nucleus. 

An  examination  of  Table  2  shows  that  in  the  B  experiment  for  a  thermal  neutron  Qux  of 
2.1  x  1017  cm^sec4  the  pumping  rate  y  creates  an  SF  photon  emission  rate  in  the  direction  of 
maximum  radiation  100  times  higher  than  that  of  the  coexisting  SE.  This  implies  that  a  thermal 
neutron  Qux  not  much  less  than  2  x  1017  is  needed  to  observe  SF  in  the  58.6  IceV  transition  of 
*®Co,  when  both  the  inhomogeneous  broadening  effect  is  small  and  the  Borrmann  effect  is  strong. 

Table  2 


Predictions  for  *°Co  Experiments 


Total  Count  Ratio 

SF/SE 

Experiment  A 

Experiment  B 

Polycrystal 

Single  crystal 
acicular  rod 

Case  1 

J  »  2.1  x  KFcm  ^sec'1 
a  =  0  b  **  0  u  =  0 

2.9 

1.9  x  10w 

Case  2 

J  *  2.1  x  lO^cnr^sec-' 
a  =  50  b  *  0  p  =  0 

104 

6.4  x  10* 

Case  3 

J  *  2.1  x  lOncm'Jsec‘l 
a  =  20  b  =  0  p  =  12  cm  -1 

8.6  x  10* 

5.5  x  104 

Case  4 

J  *  2.1  x  10,7cnr%ec-1 
a  =  0  b*0  p  =  0 

1.6  x  104 

100 

IV-  Conclusions 

The  suggested  experiment  with  *°Co  as  a  source  of  nuclear  SF  appears  to  be  feasible. 
The  thermal  neutron  pumping  flux  required  to  prepare  the  isomer  in  its  inverted  state  is  of  the 
order  of  l0,7neutrons  cm^sec'1,  which  is  within  reach  of  existing  technology3.  Analyses  similar  to 
that  presented  in  this  paper  applied  to  other  candidates,  perhaps  involving  other  pumping  meth¬ 
ods,  may  uncover  even  more  promising  nuclear  SF  sources. 
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I.  Introduction 

Superfluorescence  (SF)  (Refs.  1-3)  is  an  example  of  the  spontaneous,  cooperative 
emission  of  coherent  radiation  by  a  collection  of  identical  molecules,  atoms  or  nuclei  Intense, 
directed  pulses  of  duration  much  shorter  than  those  due  to  the  spontaneous  emission  (SE)  or 
amplified  spontaneous  emission  (ASE)  pulses  characterize  the  radiation.  The  emission  time  of  a 
cooperative  pulse  is  inversely  proportional  to  the  number  of  cooperating  emitters,  and  is  therefore 
likely  to  be  much  shorter  than  the  natural  lifetime.  Also,  the  peak  intensity  of  the  radiated  pulse  is 
proportional  to  the  square  of  the  number  of  cooperating  emitters  and  therefore  much  larger  than 
either  the  intensity  of  the  incoherent  SE,  which  is  just  proportional  to  the  number  of  emitters,  or  the 

intensity  of  ASE.  Although  the  kinetics  of  stimulated  ^radiation  has  been  discussed  in  the 

literature  (Ref.  4,  3,  6),  nuclear  SF,  which  may  be  the  most  important  process  in  a  y-ray  laser 
operation,  has  not  been  adequately  treated. 

SF  has  been  observed  for  atoms  and  molecules  (Ref.  2).  Existing  theoretical  models  (Ref. 
2)  have  successfully  explained  the  observations.  It  has  also  been  observed  experimentally  in 
atomic  and  molecular  systems  that  as  the  dephasing  rate  between  the  emitters  increases,  the  SF 
pulse  shape  changes  into  a  characteristic  amplified  spontaneous  emission  (ASE)  pulse  shape:  die 
peak  shifts  to  longer  times,  the  pulse  width  increases  and  intensity  becomes  proportional  to  the 
number  of  emitters.  However,  thus  far  no  one  has  reported  an  observation  of  nuclear  SF,  nor 
have  the  theoretical  models  used  to  analyze  the  corresponding  nuclear  phenomena  been  adequate 
for  predicting  whether  nuclear  SF  is  possible. 

In  this  paper  we  present  a  theory  of  nuclear  superfluorescence  (SF)  based  on  the  Haake- 
Reibold  model  for  the  atomic  case  (Ref.  7).  Certain  modifications  of  the  model  make  it  possible  to 
take  into  account  some  effects  that  are  more  important  in  nuclear  than  in  atomic  SF:  attenuation, 
competing  transitions,  finite  pumping  times,  and  both  homogeneous  and  inhomogeneous  line 
broadening.  Results  of  some  explicit  calculations  illustrate  die  influence  of  these  effects  on  the 
radiated  pulse.  An  extension  of  this  theory  can  be  used  to  model  ASE  and  the  transition  from  SF 
to  ASE  as  a  function  of  the  dephasing  rate  (homogeneous  and  inhomogeneous  broadening).  To 
check  the  validity  of  our  model  we  compared  calculated  pulse  shapes  with  published  experimental 
results  obtained  with  KQ:02  (Ref.  8)  in  the  SF  and  ASE  regimes  and  the  transition  region.  As  a 
final  example,  we  apply  the  theory  to  a  real  nucleus,  «Co,  created  by  thermal  neutron  irradiation 
of  39Co,  and  show  under  what  conditions  nuclear  SF  is  possible  in  this  system. 

IL  The  Haake-Reibold  Model  for  Nuclear  Superfluorescence:  Verification  and  Prediction 

Fig.  1  displays  the  Haake-Reibold  model  for  SF,  including  some  parameters  that  cannot  be 
neglected  in  the  nuclear  case.  The  diagram  on  the  left  shows  the  level  scheme  and  the  transitions 
assumed  in  this  paper.  The  Max well-B  loch  equations  on  the  right  govern  the  relation  between  die 

populations  N4,  N3,  Nj  of  die  three  levels,  the  electric  field  E1,  the  collective  polarization  R*„  and 
the  source  term  £*,  all  of  which  are  functions  of  distance  and  time.  The  real  time  and  space 
coordinates  are  normalized  to  the  superfluorescence  time  and  cavity  length,  respectively.  The 
superfluorescence  time  is  given  by 


tsF  «  8  XTo/SX^pl ; 


1  is  the  cavity  length,  X  is  the  wavelength  of  the  emitted  photon,  t  is  the  natural  radiative  lifetime, 
and  p  is  the  inversion  density  of  cooperating  nuclei. 


Figure  1.  The  Heake-Reibold  Model 

In  the  equations  t  represents  retarded  time,  the  source  term  £*  (a  stochastic  quantity) 
represents  noise  due  to  fluctuations  of  the  vacuum  the  properties  of  which  Ref.  9  derives  using  a 
quantum  electrodynanric  argument.  The  SF  model,  as  indicated  in  Fig.  1,  includes  electromagnetic 
attenuation,  natural  decay,  competing  transitions,  homogeneous  and  inhomogeneous  line 

broadening,  pumping  rate,  and  relaxation.  The  reciprocal  of  y  is  the  effective  pumping  time,  T  is 
the  natural  decay  rate  and  p  causes  attenuation  of  the  propagating  electric  field.  The  homogeneous 

line  broadening  mechanism  differs  from  that  of  die  inhomogeneous  broadening.  The  quantity  I*$, 
which  represents  dephasing  associated  with  homogeneous  line  broadening,  simply  adds  to  the 
effect  of  the  natural  decay  while  the  function  g(t)  takes  into  account  inhomogeneous  line 

broadening.  The  form  of  g(t)  depends  on  the  statistical  distribution  of  frequencies  toj  associated 
with  random  dephasing  factors  contributed  at  a  particular  location  by  photons  emitted  at  other 
positions  xj  along  the  active  region.  If  die  distribution  is  Lorentzian  g(t)  has  the  form 

g(0  *  So  e*P  (•  Tgt/2) 

go  is  the  coupling  constant  relating  the  electric  field  to  the  polarization. 

Using  a  Monte  Carlo  procedure  to  simulate  the  stochastic  source  £*,  numerical  solutions  of 
the  Maxwell-Bloch  equations  in  Fig.  1  demonstrate  both  the  attenuation  p  and  the  line  broadening 

Pe  or  1%.  Inhomogeneous  broadening  has  a  more  drastic  effect  on  the  line  shape  than 
homogeneous  broadening  for  die  same  amount  of  broadening. 

Ref.  8,  reporting  on  observations  of  combined  SF  and  amplified  spontaneous  emission 
(ASE)  produced  by  a  KQ  02  atomic  system,  gives  experimental  curves  showing  the  effect  of  line 
broadening  on  the  delay  time  of  emitted  pulses  of  both  types.  We  performed  calculations  based  on 
the  Haake-Reibold  to  compare  with  the  experimental  results.  Reasonable  values  were  chosen  for 
parameters  that  could  not  be  measured  or  otherwise  determined. 

In  the  Haake-Reibold  (Fig.  7)  model  the  stochastic  source  term  is  proportional  to  the 
pumping  rate  yand  decays  exponentially  with  the  rate  y.  This  source,  £*  generates  SF  pulses 


agreeing  with  experimental  results  in  the  SF  regime  but  not  in  the  ASE  regime.  In  the  ASE  regime 
a  stochastic  source  with  y  replaced  by  I* ,  the  spontaneous  lifetime  of  the  excited  state,  provides 
agreement  with  experiment.  Altemadvely  an  incoming  pulse  at  the  left  end  of  the  active  region 
with  an  amplitude  decaying  at  the  rate  T  can  be  used  to  model  ASE 

Figure.  2  shows  the  comparison  of  the  theoretical  calculations  with  the  experimental  data  for 
five  different  temperatures  (Ref.  8).  The  result  at  10K  which  we  assumed  to  be  due  to  pure  SF 

was  used  to  determine  the  size  of  the  active  volume  or  cooperation  number  N  and  the  SF  time  tsF- 

The  rest  of  the  calculated  results  were  obtained  using  N-109.  an  inhomogeneous  broadening  of  400 
times  the  natural  line  width  and  homogeneous  broadening  determined  by  the  respective 
temperatures  as  shown  in  the  figure. 

In  Figure  2  the  experimental  pulses  shown  are  reproduced  from  Ref.  8.  The  solid  curves 
give  the  calculated  results  obtained  with  Tsf  =  1.82  ps.  The  dashed  curves  represent  calculated 
results  with  Xsf  doubled. 


<*•  to 


Figure  2  Experimental  and  Theoretical  SF  and  ASE  pulae  altapea  from  KCIA'j 

Calculations  based  on  the  Haake-Reibold  model  used  to  predin  SF  and  ASE  agree  quite 
well  with  the  Ref.  8  results  when  reasonable  values  are  chosen  for  parameters  that  could  not  be 
measured  or  otherwise  determined. 

Unfortunately,  no  one  has  yet  reported  a  similar  experiment  involving  nuclear  rather  than 
atomic  emissions.  However,  Ref.  10  gives  a  list  of  candidates  to  replace  KC1:0"2  in  such  an 
experiment,  one  of  which  is  toCo  prepared  in  the  excited  isomeric  state  by  thermal  neutron 
pumping  of  39Co. 

Aside  from  this  pumping  mechanism,  other  advantages  of  60Co  are:  (1)  its  lifetime,  906 
sec.,  which  would  allow  long  pumping  times  for  population  inversion  and  therefore  would  involve 
achievable  thermal  neutron  fluxes;  (2)  the  fact  that  its  natural  lifetime  is,  nevertheless,  relatively 
short  and  thus  its  coupling  with  the  electromagnedc  field  relatively  strong;  (3)  its  relatively  good 
Mfiosbauer  effect;  (4)  its  low  losses  due  to  such  factors  as  internal  conversion  and  electronic 
attenuation. 
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Using  characteristic  parameters  for  60Co  and  the  thermal  neutron  cross-section  for 
preparing  *°Co  from  "Co  we  calculated  the  SF  pulse  intensity  for  different  values  of 
inhomogeneous  and  homogeneous  broadening  and  different  thermal  neutron  fluxes.  Assuming  an 
angular  shape  for  the  active  material  we  compared  the  SF  intensity  in  an  end  fire  mode  with  the 
natural  decay  component  along  the  axis.  We  found  that  for  a  thermal  neutron  flux  of  2.1  x  1017 
cnr2sec*1  the  pumping  rate  y  creates  an  SF  photon  emission  rate  in  the  direction  of  maximum 
radiation  100  times  higher  than  that  of  the  coexisting  SE.  This  implies  that  a  thermal  neutron  flux 
can  be  used  to  observe  SF  in  the  58.6  keV  transition  of  60Co,  when  both  the  inhomogeneous 
broadening  effect  is  small  and  the  Borrmann  effect  which  reduces  the  attenuation  in  certain 
directions  in  a  crystal,  is  strong.  Reactors  providing  such  thermal  neutron  fluxes  have  been 
designed  and  may  become  available  (Ref.  9, 10). 

m.  Conclusions 

The  Haake-Reibold  model  for  SF  has  been  generalized  to  treat  nuclear  SF  and  ASE.  A 
comparison  of  theoretical  results  with  experimental  pulses  obtained  from  KQ^j  gives  good 
agreement  far  atomic  SF,  ASE  and  the  transition  region.  This  agreement  with  experiment  justifies 
some  confidence  in  the  model  for  predicting  the  conditions  required  for  the  occurrence  of  nuclear 
SF. 


The  suggested  experiment  with  60Co  as  a  source  of  nuclear  SF  appears  to  be  feasible.  The 
thermal  neutron  pumping  flux  required  to  prepare  the  isomer  in  its  inverted  state  is  of  die  order  of 
10l7  neutrons  cm*2sec*1,  which  is  within  reach  of  existing  technology.  Analyses  similar  to  that 
presented  in  this  paper  applied  to  other  candidates,  perhaps  involving  other  pumping  methods,  may 
uncover  even  more  promising  nuclear  SF  sources. 
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In  a  recent  paper  (R.  Coussement  et  at.,  Phys.  Rev.  B  45,  9755  (1992)]  it  was  suggested  that  for  a 
long-lived  state  that  can  undergo  a  Mossbauer  transition,  if  the  homogeneously  broadened  width,  rather 
than  the  natural  width,  exceeds  the  solid-state-induced  inhomogeneous  width,  it  should  be  possible  to 
observe  resonant  emission  and  absorption.  On  the  other  hand,  it  is  expected  that  the  probability  of  ab¬ 
sorption  from  a  spectral  distribution  of  total  width  T  should  be  proportional  to  rr/P,  where  Tr  is  the 
radiation  width  of  the  absorbing  nucleus.  Generally  I"  is  expected  to  be  larger  than  1%  because  of  the 
possibilities  of  internal  conversion,  transitions  between  magnetic  sublevels,  etc.,  all  of  which  contribute 
to  homogeneous  broadening  and  should  reduce  the  Mossbauer  effect.  The  present  Comment  addresses 
questions  of  principle,  mathematical  assumptions,  and  other  aspects  of  the  work  of  Coussement  et  al. 
and  demonstrates  that  a  correct  statistical  treatment  recovers  the  conventional  more  stringent  require¬ 
ment  on  the  broadening. 


Coussement  et  al.1  have  argued  that  the  usual  criterion  linewidth  of  the  M5ssbauer  transition.  The  distinction 

for  the  observation  of  Mossbauer  self-absorption  is  too  between  the  natural  linewidth  r„,  the  total  linewidth  T, 

pessimistic:  'Tt  is  merely  sufficient  for  the  inhomogeneous  and  the  radiative  linewidth  Tr  is  not  made  in  Ref.  1  and 

width  (or  detuning  distribution)  to  be  small  compared  thus  the  physical  result  we  derive  here  is  notationally  for- 

with  the  width  of  the  Fourier  spectrum  of  the  fluctua-  bidden  in  Ref.  I.  In  Ref.  1  the  energy  difference  bto  is  de- 

lions,  rather  than  with  the  natural  linewidth."  However,  scribed  as  due  to  the  detuning  but,  as  will  be  discussed 

this  conclusion  is  in  error,  due  primarily  to  an  incon-  below,  a  more  general  interpretation  is  possible.  Refer- 

sistent  treatment  of  the  random  factors  in  the  wave  func-  ence  1  then  considers  an  average  rate  defined  by 

tions.  While  we  have  a  number  of  other  problems  with  pr+ar,r  .r. 

the  details  of  their  treatment,  they  are  incidental  to  our  Rr(6at  )  =  — - — —  -  ■  ^  ,  (2) 

central  point.  We  shall  therefore  accept  their  initial  as-  ^7 

sumptions  and  demonstrate  that  a  correct  statistical  where  AT  is  long  compared  to  the  fluctuation  time  scale 

treatment  recovers  the  conventional,  more  stringent  re-  of  fit)  and  gU)  but  short  compared  to  the  nuclear  life- 

quirement  on  the  broadening.  In  the  second  part  of  the  time.  Under  these  conditions  the  temporal  average  impli- 

Comment  we  demonstrate  the  distinction  between  the  Cit  in  the  ratio  of  the  above  expression  should  be 

overlap  of  spectral  lines  and  the  resonant  absorption  of  a  equivalent  to  computing  the  ensemble  average  of  PTlbto). 

spectral  line  and  show  that  contrary  to  the  contention  of  If  ,he  rale  for  intermediate  times  is  then  desired,  the  rate 

Ref.  1  complete  overlap  of  lines  does  not  guarantee  a  is  juJt  given  by  the  derivative  of  PT(bw).  We  will  there- 

strong  Mossbauer  effect.  fore  concentrate  on  a  correct  evaluation  of  the  probabili- 

According  to  Ref.  1,  Eq.  (8),  the  probability  that  a  ty  PT(Sa>).  Rewrite  PT(hto)  as 
transition  has  taken  place  is  proportional  to  the  following 
integral.  (The  amplitude  factor  is  omitted  for  simplicity.)  PTibw)= 

PT(bo)x  |/ordr/(r)g,<f)e*r*'“'‘-|:  ,  (1) 

where  is  the  energy  difference  in  the  absence  Taking  the  ensemble  average  over  the  random  fields  / 

of  the  random  fluctuation  factors  /(/)  and  g  (f)  for  the  and  g  and  assuming  them  to  be  independent  and  station- 

emitter  and  absorber,  respectively,  and  r,  is  the  natural  ary  (as  in  Ref.  1)  one  has 


JTdt'JTdtfU)fU')gU')g9tt) 


-r.n+n-iwi-n 
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/>r(8a))«  f  Tdt‘  f  rrfill7lI+C/U-i,)l 
w  0  y  0  J 

xil*lJ+c;(i-n] 

(4) 

where  /  and  g  are  (he  ensemble  averages  of  fit)  and  g  ( f ) 
[/■»£(/)],  J*£(g),  where  El  )  denotes  the  expectation 
value  over  the  random  perturbations)2  and  Cy  and  Cf  are 
the  two-point  correlation  functions.1 

Denoting  the  Fourier  transforms  of  C{  and  Cf  by  Fla) 
and  (7(ai )  and  defining  Qla)  to  be  the  convolution  of 
Fl «)  and  Cl  —  to),  Eq.  (4)  can  be  written  as 

Pr(6<j) 

_  r  da>  n,  1— IcosHSoi+ariTlc  r'T+e  ~r,T 

■*  2 s’  T2 -t-(6*)+a>)2 

(5a) 


with 

»fr(a))  =  0(a))+|/|iG(-a))+|flJF(a)) 


+  2tr6(o))|/l:lf  t! . 


(5b) 


Thus,  the  resulting  probability  is  a  weighted  sum  of  stan¬ 
dard  terms,  which  reduce  to  Lorentzians  for  large  T,  the 
first  term  persisting  even  if  the  average  values  of/(r)  and 
g(r)  are  zero.4 

Equation  (5a)  can  be  simply  evaluated  for  distributions 
of  interest.  For  example,  if  we  assume  that  the  average 
values  of  fit)  and  git )  vanish  and  the  autocorrelation 
functions  are  exponential  with  time  constant  1/A,  then 
for  large  times. 


P(6o))  =  -p 
*  $ 


r,+a 

(5a)):-Mr,l  +  A)2 


(6) 


be  inversely  proportional  to  T2. 

In  Ref.  1,  &cj  is  described  as  the  detuning  of  the  emi: 
ting  and  absorbing  states.  It  refers  to  one  pair  of  an  ab¬ 
sorber  and  emitter  with  slightly  different  resonance  ener 
gies.  The  interaction  between  the  nucleus  and  its  envi¬ 
ronment  is  different  at  different  lattice  sites  in  the  crystal 
This  provides  a  distribution  of  resonant  energies  and  i- 
the  source  of  inhomogeneous  broadening.  The  final  re¬ 
sult  for  the  probability  of  absorption  would  require  an  in¬ 
tegration  over  Sto  to  obtain  the  effect  of  inhomogeneous 
broadening.  In  this  case.  A  represents  a  random  broaden¬ 
ing  separate  from  the  fundamental  inhomogeneous 
broadening.  The  distribution  of  the  inhomogeneou 
broadening  detuning  terms  is  not  discussed  in  Ref.  . 
However,  it  is  easy  to  see  within  the  formalism  given  her. 
that  if  the  distribution  of  inhomogeneously  broadened 
states  is  of  the  same  form  as  the  correlation  function  use_ 
for  fit)  and  g(r),  then  the  result  is  of  the  same  form  u- 
given  above  with  an  increased  value  of  the  broademnr 
term: 

^toiil —  +  Amfeomo  ,  (S 

where  A^-  is  the  homogeneous  broadening  due  to  the  fluc¬ 
tuation  fit).  The  value  of  6 o>  [Eqs.  (6)  and  (7)]  after  Sud¬ 
an  integration  would  be  any  remaining  energy  different  - 
between  the  mean  energies  of  the  inhomogeneous!; 
broadened  states,  for  example,  that  given  by  any  Dopplt.' 
shift  between  emitter  and  absorber. 

On  the  other  hand,  one  can  use  the  formalism  giver 
above  to  represent  the  inhomogeneously  broadened  term 
directly.  For  a  set  of  fixed  states  with  energies  shifted  b; 
tk,  the  ensemble  average  has  the  interpretation  of  a  sur 
over  the  possible  k  states  of  the  emitter  and  absorber  anc 

k 


When  the  broadening  due  to  the  random  fluctuations,  A, 
is  small  compared  to  the  basic  linewidth,  the  unperturbed 
Lorentzian  line  shape  is  recovered.  On  the  other  hand, 
for  large  A,  the  probability  at  the  peak  5<u=0  is  inversely 
proportional  to  T,A. 

A  similar  result  is  obtained  if  a  Gaussian  autocorrela¬ 
tion  function  is  assumed.  In  this  case,  we  have1 

F(6o>)  =  vV/2— ^■■■Rc[exp(z)2erfc(z)] ,  (7a) 

r.  +/&ai 

z=-vn-  ’  (7b) 

where  A  is  the  standard  deviation  of  the  Gaussian  distri¬ 
bution  of  Ql to).  The  behavior  of  this  expression  is  more 
complicated  than  that  of  Eq.  (6)  but  it  has  the  same  limits 
for  large  and  small  A.  Similar  expressions  can  be  ob¬ 
tained  for  finite  values  of  T. 

For  both  the  exponential  and  Gaussian  distributions, 
the  price  paid  for  the  homogeneous  broadening  is  that 
the  absorption  rate  on  resonance  is  decreased  from  the 
unbroadened  case,  by  the  same  factor  as  the  increase  in 
width,  with  the  probability  inversely  proportional  to 
r„A.  In  contrast,  the  probability  given  in  Ref.  1  would 


representing  a  discrete  sum  of  individual  unbroadene- 
lines.  Note  that  this  form  of  the  correlation  function  i? 
not  necessarily  real.  Assuming  the  forest  of  lines  to  br 
essentially  dense,  using  cither  an  exponential  or  Gaussia: 
distribution  for  pltk)  and  replacing  the  sum  by  an  ir. 
tegral  reproduces  the  results  given  in  Eqs.  (6)  and  (7)  wit: 
A  representing  the  inhomogeneous  width  and  &o)  bein. 
the  difference  between  the  means  of  the  emitting  and  ab¬ 
sorbing  states  (typically,  the  Doppler  shift). 

These  additional  considerations  do  not  support  th 
Ref.  1  conclusion  that  homogeneous  broadening  due  i 
relaxation  increases  the  recoilless  self-absorption  in  tk_ 
,09Ag47  experiment,  even  though  overlap  of  lines  may  t- 
achieved.  Overlap  of  lines  is  not  a  sufficient  condition  fc 
a  strong  resonant  absorption  effect  as  the  authors  of  Re 
1  claim.  A  vanishing  effect  might  be  expected  on  pure! 
physical  grounds  in  the  limit  of  large  homogeneou 
broadening  even  as  the  overlap  between  the  emission  lir.. 
and  the  absorption  line  increases. 

An  illustrative  model  that  distinguishes  between  th 
overlap  of  spectra)  lines  and  the  probability  of  emissic 
followed  by  absorption  consists  of  two  Lorentzian  line 
with  a  common  (unperturbed)  resonance  energy  £0  and 
common  total  width  T,  which  are  displaced  from  cai 
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other  by  28.  The  overbp  of  these  two  lines  is  given  by 


A  ( r,6)  = 


dE  /,(£)  +  /  “«/£/,(£)  . 

'  *o 


(10a) 


where 

/,.,<£)« 


<r/2»r) 

(£-£0=6)J+<r/2>: 


Performing  the  integration  yields 
j4(I\8)m  1— (2/tr)tan-l(26/T)  . 


(10b) 


(11) 


in  the  limit  T— *«,  the  limit  of  large  homogeneous 
broadening,  the  area  /4(r,8)— »1  and  complete  overlap 
obtains.  This  result,  however,  does  not  imply  that  there 
is  a  strong  resonance  effect. 

The  resonance  effect  for  a  nuclear  emission  process  fol¬ 
lowed  by  an  absorption  process  is  given  by 


/nrr,r,8)“  /_"_<*£ 


/0(r/2r) 

(£-£0+6):+(r/2)J 


x  *orrr/4 

<£—£„— 8)2+(T/2)j  ’ 


(12) 


where  the  resonance  energies  are  displaced  from  each 
other  by  28.  The  first  Lorentzian  represents  the  beam 
profile  and  is  normalized  to  unity  so  that  the  beam  inten¬ 
sity  is  given  by  /0.  The  second  Lorentzian  represents  the 
scattering  or  absorption  of  a  photon  with  a  radiation 
width  for  the  process  Fr,  a  total  width  T,  and  o0  is  the 
on-resonance  cross  section  X2(2/f  + 1  )/[2ir<2/f  + 1 )( 1 
+a)j.  with  I,  (/f  1  che  excited-  (ground-)  state  spin,  the 
internal  conversion  coefficient  is  a.  and  k  is  the  nominal 
wavelength  of  the  radiation. 

Performing  the  integration  of  Eq.  (12)  yields 


£(rr.r,8)= 


A)go 

2 


1 

1-H28/D1  ' 


(13) 


Increasing  the  detuning.  8,  decreases  the  resonance  effect. 


For  a  large  detuning,  6  »  T.  and  for  minimal  homogene¬ 
ous  broadening  r  =  rr,  and  therefore  R(Vr,rr&) 
a  Fy/6: «  1.  As  the  homogeneous  broadening  increases 
the  resonance  effect  increases,  but  even  at  the  maximum, 
r-26.  is  still  small  since  R  ( rv.26,6)=s  ry/6  «  1.  In¬ 
creasing  the  homogeneous  broadening  beyond  T  —  28  de¬ 
creases  the  resonance  effect  unless  Fy  is  identical  with  F, 
i.e.,  unless  the  transition  of  interest  is  responsible  for  the 
total  width  and  Fr  is  assumed  not  to  be  constant.  This  is 
not  a  realistic  assumption,  however,  since  in  nuclear  sys¬ 
tems  of  interest  the  total  width  is  usually  dominated  by 
internal  conversion  and  other  processes  while  the  radia¬ 
tive  width  is  only  a  small  part.  Not  making  a  distinction 
between  Tr  and  F  in  Eq.  (13)  has  the  same  effect  as  using 
the  overlap  in  Eq.  (11)  to  calculate  this  resonance  effect. 
Alternatively,  introducing  Tr.  for  T  in  the  numerator  of 
(10b)  reduces  the  overbp  in  Eq.  (11)  by  the  factor  rr/r. 
Since  we  can  hardly  change  the  radiative  width  Tr  of  a 
nuclear  transition  but  we  can  easily  change  the  total  or 
homogeneous  width,  it  is  important  to  keep  these  distinc¬ 
tions  dear  in  the  analysis. 

To  model  the  inhomogeneous  broadening  the  integral 
over  appropriate  distribution  of  the  detuning  parameter  8 
in  Eq.  (13)  must  be  taken.  In  another  paper5  the  result  is 
shown  to  be  the  error  function  of  Eq.  (7a)  multiplied  by 
the  appropriate  factor.  When  magnetic  sublevels  may  re¬ 
lax  among  themselves,  increasing  homogeneous  broaden¬ 
ing  increases  the  relaxation  rate  and  F  but  does  not  alter 
Tr  In  the  limit  of  infinite  homogeneous  broadening,  the 
tout  width  T—  cc ,  and  although  the  last  factor  on  the 
right  in  Eq.  (13)  tends  to  unity,  the  ratio  ry/r— 0.  This 
result  remains  valid  when  the  present  treatment  is  gen¬ 
eralized. 

In  another  paper,3  we  present  a  general  treatment  of 
the  emission  and  absorption  of  photons  by  nuclei  with  in- 
homogeneously  broadened  linewidths.  This  treatment, 
based  on  the  results  of  Boyle  and  Hall"  and  Heitler7  is 
more  appropriate  for  the  analysis  of  Mossbauer  experi¬ 
ments  with  isomeric  nuclei. 


*R.  Coussement.  G.  S’heeren,  M.  Van  Den  Bcrgh,  and  P.  Bool- 
chand.  Phys.  Rev.  B  4S.  975S  (1992). 

2Note  that  Ref.  1  has  an  inconsistency  in  the  treatment  of  El/) 
and  Elg).  It  is  implied  that  the  average  of  /  and  that  the 
average  of  the  magnitude-squared  of  fit)  are  both  unity 
{£(/)»■  1  and  £(//’)«!].  This  cannot  hold  unless  /  is 
identically  one.  If  one  were  to  assume  as  in  Blume  and  Tjon 
{Phys.  Rev.  J65,  446  (1968)]  that  fit)  and  git)  were  random 
phase  functions,  i.e.,  that  /(r)Kexp[i4U)]>  the  latter  condi¬ 
tion  would  hold  identically  and  one  might  expect  £(/)“0. 
For  the  purposes  of  this  comment  the  precise  value  of  £  If)  is 
unimportant. 

3One  expects  the  correlation  functions  to  be  real  for  random 


homogeneous  broadening  but  there  is  no  loss  in  assuming 
them  to  be  complex;  the  difference  expressed  in  Eq.  (4)  corre¬ 
sponds  to  the  difference  between  emitting  and  absorbing  sites. 

4If  the  correlation  functions  C/  and  C ,  are  real.  F,  C,  and  Q  are 
even.  In  any  event,  F,  C,  and  Q  are  always  real.  The  use  of 
Gl  —  u)  in  contrast  to  £(o)  reflects  in  Ref.  3  the  effect  at  em¬ 
itting  and  absorbing  sites. 

3B.  Balko,  L.  Cohen,  D.  A.  Sparrow,  and  J.  F.  Nicoll  (unpub¬ 
lished).  See  also  B.  Balko,  G.  Heriing,  I.  W.  Kay,  and  J. 
Nicoll  (unpublished). 

‘A.  J.  F.  Boyle  and  H.  E.  Hall,  Rep.  Prog.  Phys.  25,  44  (1962). 

7W.  Heitler,  The  Quantum  Theory  of  Radiation,  3rd  ed.  (Ox¬ 
ford  University  Press,  London,  1954). 
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The  physical  assumptions  made  to  calculate  resonant  absorption  for  the  nuclear  transitions  from 
long-lived  isomeric  (r>  1  ms)  states  by  Coussement  et  al.  [Phys.  Rev.  B  4S,  9733  (1992)]  differ  qualita¬ 
tively  from  those  made  to  calculate  nuclear  resonant  absorption  for  transitions  from  the  conventional 
short-lived  (r  <0. 1  ms)  states.  It  is  for  this  reason  that  one  cannot  mathematically  deduce  the  resonant 
fraction  for  the  case  of  the  long-lived  states  by  mere  extrapolation  of  the  formulas  from  the  case  of 
short-lived  ones. 


In  our  paper1  we  have  chimed  that  the  Mossbauer 
effect  in  long-lived  nuclear  states  (r>  I  s)  will  be  possible 
if  the  inhomogeneous  width  (T^)  is  less  that  the  homo¬ 
geneous  or  relaxation  width  ( T,,  rather  than  the  natural 
width  (r„ ).  Our  claim  in  this  context  concerns 
specifically  (a)  long-lived  states  and  (b)  with  both  the 
source  and  absorber  matrices  homogeneously  broadened 
by  fluctuations.  In  Mossbauer  spectroscopy  with  conven¬ 
tional  short-lived  nuclear  states  (r<l  /is),  it  is  well 
known  that  to  observe  the  resonance  effect,  the  inhomo¬ 
geneous  width  rinh  should  be  less  than  I*,,  representing  a 
more  stringent  requirement,  which  we  clearly  accept.  We 
believe  that  in  their  Comment,  Balko  et  al.2  argue  that 
our  claim  of  the  less  stringent  requirement,  viz., 
r,rt  <  r„  to  observe  the  Mossbauer  effect  in  long-lived 
states  is  actually  in  error  due  to  an  inconsistent  treatment 
of  random  factors  in  the  wave  function.  Balko  et  al.1 
claim  that  if  the  algebra  is  done  right,  one  recovers  the 
more  stringent  requirement,  viz.,  <  T,  to  observe  the 
Mossbauer  effect  in  long-lived  isotopes  as  well.  We 
disagree  with  this  conclusion.  We  would  like  to  address 
the  issue  in  two  parts,  at  first  conceptually  and  then 
mathematically. 

Conceptually,  the  requirement  to  observe  the 
Mossbauer  effect  in  long-lived  states  differs  qualitatively 
from  the  case  of  short-lived  ones,  because  of  the  widely 
different  temporal  behavior  of  the  resonant  absorption 
process  in  the  two  cases.  In  the  long-lived  case  (r  — 
seconds),  as  we  indicated  in  our  paper,  the  interaction  of 
the  electromagnetic  (EM)  field  with  nuclei  (ff/r  =* 
10~,7eV)  becomes  much  weaker  than  the  interaction  of 
nuclei  with  its  environment  (fi/T,=*10_M  eV)  in  a  crys¬ 
talline  solid.  This  implies  that  during  the  long  time  that 
nuclear  resonant  absorption  occurs,  nuclei  in  a  crystalline 
solid  will  have  ample  opportunity  to  exchange  energy 
with  their  environment,  and  one  expects  the  detuning 
width  or  resonance  width  to  be  homogeneously 
broadened  by  a  characteristic  relaxation  time  r,  (of  order 
of  few  ms)  of  the  ensemble.  The  absorption  effect  will  be 
observed  if  ruUl<r,. 
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In  the  case  of  short-lived  nuclear  states  (r  <  /is),  the  :: 
quirements  on  times  and  therefore  energies,  are  complete 
ly  reversed.  The  interaction  energy  of  the  EM  field  wr..- 
nuclei  (a  10-M>  eV)  far  exceeds  the  interaction  energy  : 
nuclei  with  their  environment  (=r  10-u  eV)  in  a  crvsi- 
line  solid.  Generally  speaking,  the  process  of  resonar, 
absorption  can  now  be  viewed  to  be  instantaneous  in  reL 
tion  to  the  relaxation  times  in  the  absorbing  medium 
This  has  the  important  consequence  that  in  resonant  ar 
sorption,  the  EM  field  encounters  a  snapshot  of  all  nucic 
present  in  a  variety  of  local  fields  in  the  crystal.  The  re- 
onance  effect  is  crucially  sensitive  to  inhomogeneou 
broadening  in  Che  lattice.  This  requires  the  me 
stringent  condition,  viz.,  ruA<r,.  We  are  thus  led  ;. 
believe  that  the  resonant  absorption  process  from  Ion: 
lived  nuclear  states  differs  qualitatively  from  short-live- 
ones.  To  observe  these  processes  it  is  thus  reasonabi. 
that  the  same  stringent  requirement,  i.e.,  rmh<r.,  ma 
not  be  necessary  in  both  cases  as  Balko  et  al.  hav? 
claimed7  in  their  Comment. 

Mathematically,  crucial  differences  come  in.  First  c 
all,  one  needs  to  bring  in  fluctuations  [fit)  and  g</))  in  - 
natural  way  to  describe  a  condition  of  equilibrium  for  thr 
case  of  long-lived  nuclear  states.  In  our  paper,1  we  chosr 
to  define  the  functions  fit)  and  g  ( t )  as  generally  as  poss: 
ble.  It  would  be  incorrect  to  talk  about  frequency  an. 
phase  modulation  of  the  wave  function  independent; 
when  clearly  one  type  of  modulation  induces  the  other 
In  our  description  of  the  fluctuations,  it  would  then  be  in¬ 
correct  to  assume  that  a  random  phase  diffusion  occur, 
because  this  would  imply  a  decaying  nuclear  state,  wher. 
actually  what  is  being  implied  is  that  changes  in  the  loca 
chemical  environment  occur,  inducing  changes  on  th. 
fields  sensed  by  the  nuclei. 

Of  major  importance  is  the  fact  that  Balko  et  al.  :  in¬ 
troduce  correlation  functions.  Such  correlation  function 
can  be  used  to  calculate  ensemble  averages  for  the  case  o: 
short-lived  states.  Such  a  procedure  for  the  case  of  long- 
lived  nuclear  states  poses  difficulties,  however.  It  is  clear 
that  the  ensemble  average  of  the  transition  probabilit: 
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will  require  calculating  an  average  over  some  integrals. 
In  their  Comment,  Balko  et  al.  proceed  further  by  as¬ 
suming  that  the  average  over  the  integral  is  equal  to  the 
integral  over  the  average  as  is  normally  done  for  short¬ 
lived  isotopes.  This  automatically  leads  to  the  definition 
of  correlation  functions.  Such  a  procedure  can  only  be 
limited  to  the  case  of  short-lived  nuclear  states  as  we  dis¬ 
cuss  next. 


Specifically.  Balko  et  al.  give  an  expression  for  the 
probability  in  their  Eq.  (3)  which  reads 

fjdt'jjdt flt)fmlt')gU')g*it) 

- r  i&4*r  —  /■) 

(1) 

In  obtaining  the  ensemble  average  of  the  probability 


fjdt'  fjdt  f{t>flt’)g(t’)g*{t)e~r‘U~'  '  'j 

in  general,  one  cannot  equate  the  average  over  the  integrals  to  the  integral  over  the  average. 
<f,r(&&»)>  =  fjdt' /o7d»^/U)/•U')g(r')g•(r)je_^",,^,  . 

I - 


But  exceptions  occur,  and  the  case  of  short-lived  states  is 
one  of  them.  Indeed  for  short-lived  states,  the  exponen¬ 
tials  t  *  and  e  *  are  functions  which  decay  much 
faster  in  time  than  the  fluctuations  fit)  and  git).  There¬ 
fore.  contributions  to  the  integrals  in  (2)  over  r  and  over 
t'  are  always  limited  to  very  short  time  intervals.  During 
such  short  time  intervals  the  fluctuation  functions  will 
not  change  appreciably  and  each  absorption  event  will 
capture  a  different  moment  and  a  different  possible  value 
of  these  functions.  In  such  a  case  one  can  reduce  the  en¬ 
semble  averaged  probability  to  an  integration  over  the  en¬ 
semble  average  of  the  fluctuating  factors  by  introducing 
correlation  functions,  as  done  by  Balko  et  a/.2  and  oth¬ 
ers. 

For  long-lived  isotopes,  however,  e  '  and  e  '  de¬ 
cay  so  slowly  that  during  the  long  emission  or  absorption 
event,  fit)  and  g(r)  can  fluctuate  over  the  full  range  of 
parameter  values  that  characterize  such  fluctuation  func¬ 
tions.  The  procedure,  followed  for  short-lived  States,  viz., 
sampling  the  fluctuating  factor  at  the  exact  transition 
moment,  can  no  longer  be  followed  for  long-lived  states 
because  each  transition  now  involves  a  long  temporal 
process  compared  to  the  fluctuation  rates.  For  long-lived 
states,  the  integrations  in  (2)  will  intrinsically  result  in  a 
time  averaging  and  each  nuclear  event  will  contain  the 
same  time  averaged  information  over  these  fluctuations. 


Because  we  can  assume  that  time  averaging  is  equal  to  an 
ensemble  averaging  at  a  given  moment,  there  is  no  fur¬ 
ther  need  to  calculate  the  ensemble  average  for  the  long- 
lived  case. 

Finally,  if  one  considers  the  situation  when  a 
Mossbauer  resonant  absorption  experiment  is  performed 
using  a  nearly  monochromatic  source  (i.e.  no  fluctuations 
in  the  source  matrix)  to  scan  an  absorber  in  which  fluc¬ 
tuation  broadening  occurs,  then  we  do  expect  the  max¬ 
imum  in  the  resonant  effect  to  decrease.  This  can  be  seen 
in  Eq.  (19)  of  Ref.  1,  where  if  we  replace  the  fATio)—b>0) 
function  by  a  b*ria>— uj  function,  the  overlap  between 
the  source  and  absorber  line  shape  decreases.  The  re¬ 
verse  situation,  i.e.,  when  the  source  matrix  is  fluctuation 
broadened  and  the  absorber  matrix  nearly  monochromat¬ 
ic,  will  likewise  also  give  rise  to  a  reduced  maximum  in 
resonant  effect.  However,  if  both  the  source  and  absorber 
matrices  are  fluctuation  broadened  in  an  uncorrelated 
manner,  then  our  claim  in  Ref.  1  is  that  the  maximum  in 
the  resonant  effect  will  be  restored. 

In  summary,  we  believe  that  the  Comment  of  Balko 
et  al.  pertains  to  the  usual  case  of  the  Mossbauer  effect 
from  short-lived  nuclear  states  ( r  £  lfis).  However  these 
calculations,  in  our  opinion,  cannot  be  extended  to  the 
case  of  the  Mossbauer  effect  from  long-lived  nuclear 
states  for  reasons  discussed  above  and  in  our  paper.1 


•R.  Coussement,  G.  S’heeren,  M.  Van  Den  Bergh,  and  P.  Bool- 
chand.  Phyv  Rev.  B  45, 97SS  (1992). 

JB.  Balko,  I.  W.  Kay.  J.  F.  Nicoll,  J.  D.  Silk,  D.  A.  Sparrow, 


and  G.  H.  Heriing,  preceding  paper,  Phys.  Rev.  B  48,  16  139 
(1993). 
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Title:  Reply  to  “Reply  to  ‘Comment  on  Nuclear  resonant  absorption  in  long-lived 
isomeric  transitions’"  by  R.  Coussement,  G.  S’heeren,  M.Van  en  Bergh,  and  F. 
Boolchand 


In  their  reply1  to  our  comment2  on  their  paper3,  the  objections  of  the  authors  of 
References  1  and  3  were  based  on  a  conceptual  argument  and  a  mathematical  concern.  In 
simple  terms,  the  essence  of  the  disagreement  is  revealed  by  the  answers  to  two  questions. 

(1)  Conceptual  question:  Do  long-lived  and  short-lived  nuclear  states  require 
different  treatment  for  describing  the  emission  and  resonance  absorption  processes?  We 
say  no  for  the  problem  at  hand. 

(2)  Mathematical  question:  Can  the  average  over  the  integral  be  equated  to  the 
integral  over  the  average?  We  say  yes,  except  for  an  integrand  that  is  a  mathematically 
contrived,  pathological  stochastic  process  of  no  interest  here. 

In  what  follows  we  provide  a  more  detailed  explanation  of  our  position. 

(1)  Conceptual  question 

It  has  been  known  for  some  time  that  the  nuclear  lifetime  affects  individual  level 
widths,  but  as  far  as  fluctuations  are  concerned,  it  is  an  irrelevant  parameter  except  when 
dealing  with  nonstationary  stochastic  processes.  [4.  M.  Blume,  Hyperfine  Structure  and 
Nuclear  Radiations,  eds.  E.  Matthias  and  D.A.  Shirley,  (North  Holland,  Amsterdam, 
1968),  p.  91 1].  Blume's  theoretical  treatment  is  completely  general,  as  far  as  the  relevant 
parameters  are  concerned.  Both  short  and  long  lived  Mbssbauer  isotopes  are  covered  by 
the  theory.  Consequendy,  the  conceptual  criticism  of  Reference  1  is  irrelevant  to  the  issue. 
If  there  is  any  substance  to  the  claim  in  References  1  and  3,  i.e.,  that  long-lived  and  short¬ 
lived  states  require  different  treatment  for  describing  die  emission  and  resonance  absorption 
processes  in  the  presence  of  fluctuating  fields,  it  has  to  be  found  in  the  mathematical 
analysis. 


(2)  Mathematical  question: 


The  mathematical  cridcism,  which  is  the  essence  of  the  Reference  1  argument,  is 
both  confused  and  wrong.  It  is  true  that  the  time  evolution  of  the  wave  function  is 
determined  by  the  usual  quantum  mechanical  expression 


V(0  =exp 


dt'  H(t') 


V(0). 


34 


It  is  the  time  dependent  H(t')  that  contains  the  stationary  random  process  implicit  in 
Reference  3.  Although  that  paper  characterizes  the  wave  function  as  non-stationary 
because  of  the  factor  e*rnt/2t  it  assumes  that  the  stochastic  process  represented  by  the  factor 
f(t)  is,  in  fact,  stationary  in  connection  with  equation  (2)  of  Reference  3.  The 
representation  of  this  time  dependent  operator  as  a  factor  f(t)  in  the  wave  function  defined 
by  equation  (1)  of  the  paper,  i.e.,  in 

-r„i/2  -»<y 

v(x,  t)  -  v(x,  0) /(t)e  "  e 

is  an  approximation  whose  validity  must  be  determined  in  each  problem.  In  general,  f(t) 
need  not  be  an  example  of  a  stationary  random  process  since  it  may  depend  on  the 
collective  state  of  all  of  the  nuclei.  However,  the  assumptions  of  Reference  3  (cf.  equation 
(2)  and  the  surrounding  discussion)  are  equivalent  to  at  least  a  weak  ft, .in  of  stationarity. 
Given  stationarity  of  f(t),  the  approach  used  in  Reference  2  is  correct  unless  the  function 
f(t)  is  pathological  in  nature. 

Reference  1  ’s  assertion  that  “in  general,  one  cannot  equate  the  average  over  the 
integral  to  the  integral  over  the  average”  is  not  applicable  to  the  present  case.  On  the 
contrary,  it  is  well  known  in  the  theory  of  stochastic  processes  that  under  very  general 
conditions: 

If  A  is  a  measurable  parameter  set  and  if  the  average  E{x  t  (to))  of  x  t  (co)  exists  for  teA, 
and 


<oo. 


then 


J  xt(Q))dt  | »  j  E  {jcf(ct))  ytx. 


[5.  J.L.  Doob,  Stochastic  Processes,  (John  Wiley  &  Sons,  New  York  1953),  pp.  62  II]. 
In  other  words,  for  absolutely  convergent  iterated  integrals,  the  result  is  independent  of  the 
order  of  integration.  In  less  formal  terms,  the  average  of  an  integral  can  be  taken  over  the 
integrand  before  the  integration  takes  place,  except  for  pathological  cases.  There  is  no 
evidence  that  suggests  that  the  present  example  is  such  a  case;  however,  if  it  were,  the 
method  proposed  by  R.  Coussement,  et  al.  would  similarly  fail. 

Finally,  the  argument  presented  in  the  second  to  last  paragraph  of  Reference  1  again 
ignores  the  distinction  between  overlap  of  lines  and  the  resonant  effect  and  is  refuted  by 
equation  13  in  Reference  2  and  the  subsequent  discussion. 
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